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INTRODUCTION 


The  B-l  aircraft,  developed  by  the  B-l  Division  of  Rockwell 
International,  utilizes  an  Electrical  Multiplex  System  (EMUX)  for: 

- processing  and  transfer  of  serial-digital  and  discrete  data 
throughout  the  aircraft 

- power  control  of  most  electrical  loads 

- automatic  electrical  load  management 

There  are  two  independent  EMUX  systems  in  the  aircraft  for  redundancy  - 
each  having  its  own  configuration  and  software.  EMUX  is  programmed 
using  an  IBM  370  hosted  compiler.  Input  to  the  compiler  consists  of 
logic  equations,  assembly  language  instructions,  and  EMUX  system  con- 
figuration data.  EMUX  can  be  programmed  to  output  over  two  thousand 
signals  as  functions  of  over  five  thousand  input  signals.  This 
magnitude  of  equations  and  signals  complicates  the  understanding  and 
analysis  of  the  EMUX  software. 

In  order  to  improve  the  Air  Force  EMUX  software  analysis  capability, 
several  in-house  computer  routines  were  developed  from  July  1975  to 
September  1977.  All  of  these  routines  use  copies  of  the  same  data  files 
created  and  maintained  by  the  B-l  Division  for  use  in  generating  the  EMUX 
flight  software.  All  of  the  routines  were  developed  on  Aeronautical 
Systems  Division's  CDC  6600  Computer  System  using  FORTRAN  EXTENDED. 

One  aSD  Technical  Report  and  four  Avionics  Directorate  (ENA)  Engin- 
eering Reports  have  been  written  to  document  the  various  computer 
programs : 

"B-l  EMUX  Data  Tape  Conversion  Software" 

ASD-ENA-77-20 

"B-l  EMUX  Usage  Analysis  Software" 

ASD- ENA-78-2 

"B-l  EMUX  Signal-Signal  and  Signal-Box  Relationship 
Analysis  Software" 

ASD- ENA-7 8-3 

"B-l  EMUX  Logic  Equation  Regeneration  Software" 

ASD-ENA-78-4 

"B-l  EMUX  Automated  Logic  Diagrammer  Demonstration  Program" 
ASD-TR-78-17 
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This  report  discusses  the  FORTRAN  software  written  for  use  on 
the  CDC  6600  to  demonstrate  the  feasibility  of  automatically  creating 
logic  diagrams  of  boolean  equations  directly  from  the  contractor 
maintained  computer  data  files.  Diagrams,  while  being  logically 
equivalent  to  equations,  are  pictorial  representations  of  the  equa- 
tions and,  as  such,  facilitate  understanding.  However,  logic  diagrams 
are  quite  time-consuming  to  draw  manually;  estimates  upwards  of  two  or 
more  man-years  for  a single  version  of  EMUX  software  have  been  given. 
Thus,  a computerized  diagramming  technique  which  did  not  require 
special  input  data  could  provide  coat  and  time  savings  if  the  automated 
diagrams  were  of  equal  quality  (accuracy,  readability,  layout,  etc.) 
as  manual  diagrams.  Program  BOLD  (B-l  (One)  Logic  Diagrammer)  is  the 
result  of  this  feasibility  demonstration. 
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PROCRAM  PKSICN  CROUNURDLES 


Program  BOLD  was  writ  Con  in  FORTRAN  EXTENDED  for  use  on 
Aeronaut teal  Systems  Division's  CDC  t>600  computer  system  using  t he 
to l lowing  ground  rules: 

- the  layout  of  the  diagrams,  expec tally  the  positioning  of  the 
logic  operator  symbols,  must  be  of  equal  or  better  quality  when 
compared  to  manually  drawn  diagrams. 

- any  valid  logic  equation  for  KMl'X  must  be  diagranuuable,  including 
equations  with  time  delays. 

- the  diagrams  must  accurately  reflect  the  boolean  equation  logic. 

- preparation  of  special  input  data  must  be  minimised. 

- program  sice,  run  time,  and  resultant  cost  must  be  less  than 
manual  costs. 

All  of  these  ground  rules  have  been  satisfied  by  BOLD  as  will  be  discussed 
in  the  following  paragraphs. 

The  proper  layout  of  the  diagram  is  one  of  the  most  important  measures 
of  the  success  or  failure  of  any  computerised  diagramming  technique. 

Proper  layout  must  consider  the  location  of  each  of  the  logic  operator 
svmhols  (.gates)  and  the  interconnection  of  the  gates.  A major  portion  of 
the  logic  in  BOLD  is  involved  in  diagram  layout.  Listed  below  are  the 
diagram  layout  rules  which  were  established: 

- flow  of  logic  will  be  from  left  to  right. 

- no  overlap  of  gates. 

- gates  will  be  arranged  in  columns  (for  readability  and  simplicity) 
and  will  be  assigned  to  columns  from  right  to  left  to  minimise 
length  of  Interconnect  lines. 

- Interconnect  lines  will  be  straight  whenever  possible  and  will 
have  at  most  two  turns. 

- no  feedback  or  latch  interconnect  lines  will  be  used  but  will  be 
labeled  as  operator  inputs  (this  is  not  inconsistent  with  KMDX 
operation  but  is  probably  not  how  the  diagram  would  be  done  manually). 

- crossover  of  Interconnect  lines  will  be  minimised. 

- time  delays  will  be  treated  as  two  input  operators.  One  Input  is 
the  quantity  to  he  delayed  and  the  other  input  is  the  time  delay 
durat ion. 
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Figures  1 and  2 are  sample  B01 D printouts  of  two  different  equation 
diagrams.  As  can  be  seen,  the  diagrams  are  "drawn"  by  the  standard 
line  printer  rather  than  by  a continuous  line  plotter  (e.e.,  a CALCOMP 
plotter).  This  type  of  output  was  selected  primarily  because  BOLD  was 
Intended  to  be  a feasibility  demonstration  program  rather  than  a final 
product  ion-oriented  program.  As  such,  diagram  output  via  line  printer 
is  simpler  to  obtain  in  terms  of  software  logic  complexity  and  program 
turnaround  time.  A logical  modification  to  improve  diagram  readability 
to  BOLD  for  production  vise  would  be  the  conversion  to  CAl.COM!’- typo 
diagrams.  While  this  would  require  the  main  program  to  be  changed 
(approximately  25  to  30X  rewrite),  the  methodology  would  be  the  same 
and  no  other  subroutine  changes  should  be  required.  At  the  same  time, 
logic  for  feedback  lines  could  be  Included  if  deemed  desirable. 

All  equations  diagrammed  to  date  have  accurately  reflected  the  input 
equat ion  logic.  The  only  problems  encountered  but  not  solved  were  due  to 
either  improper  input  equation  syntax  or  equation  sire  in  excess  of  BOLD 
capacities.  Improper  equation  syntax  (e.g.,  missing  operator,  missing 
parenthesis)  would  also  cause  problems  for  the  contractor's  EMUX  compiler 
and  can  only  be  resolved  by  correcting  the  equation.  Many  data  arrays 
are  used  in  BOLD  for  storage  of  various  types  of  equation  information. 

The  sires  of  these  arrays  are  related  to  each  other  to  some  extent  but 
are  somewhat  arbitrarily  established.  There  are  no  known  theoretical 
limits  to  the  equation  sire  which  can  be  processed;  however,  practical 
limits  of  computer  memory  available  and  desired  turnaround  time  may 
prevent  some  "extreme"  equations  from  being  diagrammed . Present  array 
sires  have  rejected  less  than  1/2  of  IX  of  all  equations  input.  Array 
sire  limitations  would  probably  be  improved  by  the  conversion  to  CALCOMP- 
type  diagrams.  The  present  program  is  limited  to  a maximum  of  8 pages  of 
printout  per  equat ion  which,  when  cut  and  taped  together,  will  show  the 
entire  diagram.  Limits  on  particular  equation  size  parameters  (e.g., 
number  of  parenthesis  pairs,  number  of  operators,  number  of  operands  per 
operator)  are  documented  in  the  program  listings. 

BOLD  uses  the  same  data  files  maintained  by  the  contractor  for  EMUX 
flight  software  programming.  The  only  additional  input  required  is  the 
user's  selection  of  the  equation  set  to  be  diagrammed.  The  user  has  the 
flexibility  to  select  individual  equations  or  entire  subsystems  for 
diagramming. 

While  final  conclusions  on  relative  costs  (automated  diagrams  versus 
manual  diagrams)  cannot  be  made  until  a production  version  of  BOLD  is 
created,  a preliminary  cost  analysis  resulted  significantly  in  favor  of 
automated  diagrams. 


PRINTOUT  DESCRIPTION 


While  the  BOLD  printout.  Figures  1 and  2,  seems  at  first  cryptic. 

It  is  understandable  if  examined  in  parts.  The  first  line  identifies 
the  OCNEK  signal  designator  (an  aircraft  standard  signal  identification 
scheme)  for  the  equation  diagrammed,  the  aircraft  effectlvity,  and  the 
section  of  EMUX  (left  or  right)  involved.  (For  Figure  1,  the  equation 
output  signal  designator  is  2821-008  and  the  equation  is  used  in  the 
left  section  of  EMUX  on  Aircraft  3.)  The  next  three  lines  on  the  left 
identify  the  aircraft  system,  subsystem,  and  subsystem  to  which  the 
signal  is  assigned.  (For  Figure  1,  sub-subsystem  2821  is  the  Internal 
Transfer  portion  of  the  Fuel  Distribution  System.)  The  line  to  the  right 
of  the  system  labels  is  an  abbreviated  signal  description.  (For  Figure  1, 
the  signal  is  a power  control  signal  for  the  tank  1 transfer  pump  number 
2821PP1.)  The  next  set  of  lines  is  the  actual  boolean  equation  to  be 
diagrammed  as  input  to  BOLD.  Following  the  equation  is  a table  of  all 
signals  that  are  used  in  the  equation.  Included  with  each  signal 
designator  is  the  signal  description.  The  final  part  of  the  printout 
is  the  diagram.  Due  to  line  printer  limitations,  the  diagram  is  not  as 
readable  as  possible.  Each  logic  operator  (gate)  is  indicated  with  a box 
of  asterisks  with  a letter  inside.  AND  gates  use  the  letter  A,  OR  gates 
use  the  letter  0,  EXCLUSIVE  OR  gates  use  the  letter  X,  Type  1 time  delays 
use  a dollar  sign  ($),  and  Type  2 time  delays  use  a question  mark  (?). 
Gates  with  more  than  two  inputs  have  the  left  column  of  asterisks  extended 
as  necessary  to  accomodate  all  inputs.  NOT'ed  gate  inputs  are  indicated 
with  the  letter  0 in  place  of  an  asterisk  in  the  gate's  leftmost  column. 
Gate  interconnecting  lines  use  dashes,  periods,  and  letter  l's  for  their 
horizontal  parts,  corners,  and  vertical  parts  respectively.  The  reada- 
bility of  these  line  printer  diagrams  can  be  improved  significantly 
merely  by  drawing  over  the  connecting  lines  manually  - as  is  done  in 
Figure  2. 
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Figure 
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PROGRAM  OPERATION 


The  process  used  by  BOLD,  as  shown  in  Figure  3,  starts  with  reading 
the  next  equation  to  be  diagrammed  into  array  EQN.  In  order  to  simplify 
processing  throughout  the  rest  of  the  program,  the  equation  is  converted 
into  an  internal  coded  integer  format  and  is  stored  in  array  EQN1. 

Array  VAR  is  loaded  with  signal  names  appearing  in  the  equation.  The 
integer  code  is  shown  below: 


VALUES 


REPRESENT 


1 to  1000 


index  of  Ignal  name  stored  in  VAR 


1001  to  2000 

value  -1000  is  the  index  of  i 
array  OPR 

2001 

( 

left  parenthesis 

2002 

) 

right  parenthesis 

2003 

* 

AND  operator 

200 A 

0 

EXCLUSIVE  OR  operator 

2005 

+ 

OR  operator 

2006 

- 

EQUAL  operator 

2007 

f 

NOT  operator 

2010 

$ 

TYPE  1 TIME  DELAY  operator 

2011 

? 

TYPE  2 TIME  DELAY  operator 

9999 

replaced  character 

-2006 

NOT  EQUAL  operator 

-2000  to  -1001 


NOT’d  operator,  absolute  value  -1000  is  the 
index  of  an  operator  in  OPR 


-1000  to  -1 


NOT'd  operator  input,  absolute  value  is  the 
index  of  a signal  name  stored  in  VAR 


Due  to  parentheses  and  operator  precedence,  the  equation  cannot  be 
simply  scanned  from  left  to  right  to  determine  the  correct  diagram. 


k 


Therefore,  BOLD  first  scans  EQN1  for  parenthesis  pairs  and  stores  infor- 
mation on  each  pair  in  array  PAR.  Once  the  parentheses  have  been  found, 
then  the  logic  within  each  pair  is  sent  to  subroutine  PARSE  for  analysis  - 
the  innermost  pair  being  done  first.  As  each  pair  is  done,  the  logic 
and  parentheses  are  replaced  in  EQN1  with  the  resultant  output  operator's 
index  in  OPR  and  with  9999's.  After  all  parentheses  have  been  removed, 
the  entire  equation  is  sent  to  PARSE  for  a final  analysis. 

PARSE  semantically  analyzes  each  portion  of  the  equation,  taking 
irto  account  operator  precedence,  by  scanning  left  to  right  five  times 
(once  for  each  operator  type) . The  operator  precedence  used  by  PARSE  is 
Type  1 Time  Delays,  Type  2 Time  Delays,  AND  gates,  EXCLUSIVE  OR  gates, 
and  OR  gates.  During  each  of  the  five  passes,  the  appropriate  operator 
is  tested  for  and,  if  found,  is  stored  in  array  OPR  along  with  its 
associated  inputs.  Thus,  when  PARSE  is  finished  with  the  equation,  all 
of  the  semantic  equation  information  is  stored  in  array  OPR  as  shown 
below  for  operator  i: 

OPR  (i,  1) 


OPR  (i,  2) 
OTR  (i,  3) 


OPR(i,  4) 


OPR  (i,  5) 
through 
OPR  (i,  34) 


operator  type  code  (2003,  2004,  2005,  2010, 
2011) 

number  of  Inputs  to  this  operator 

x of  x,  y coordinates  of  location  of  operator 
on  output  page 

y of  x,  y coordinates  of  location  of  operator 
on  output  page 

operator  input  indices.  Negative  values  indi- 
cate primed  inputs.  Absolute  values  from  1 to 
1000  are  signal  name  indices  pointing  to  array 
VAR.  Absolute  values  from  1001  to  2000  are  in- 
puts from  other  operator  outputs.  These  values 
-1000  point  to  array  OPR. 


The  present  size  of  array  OPR  limits  any  one  operator  to  30  inputs  maxi- 
mum. It  should  be  noted  that  OPR  (i,  3)  and  OPR  (i,  4),  operator  location 
Information,  are  not  defined  by  BOLD  at  this  point  in  the  equation 
processing. 

The  only  remaining  task  before  printing  out  the  diagram  is  the 
determination  of  the  location  of  each  of  the  operators  on  the  output 
page.  This  portion  of  the  logic  within  BOLD  was  designed  to: 

- minimize  crossover  of  interconnecting  lines. 

- minimize  the  number  of  operators  which  could  not  be  connected 
due  to  the  location  of  other  operators 
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- maximise  the  uae  of  straight  line  connections 


- make  all  connections  so  that  the  logic  flow  Is  from  left  to 
right 

For  both  readability  and  ease  of  programming  it  was  decided  that  each 
operator  would  be  assigned  to  one  of  several  available  columns.  Once 
an  operator  Is  assigned  to  a column  then  its  vertical  position  In  the 
column  Is  determined  based  on  providing  the  straightest  connection 
from  Its  output  to  the  next  operator’s  input.  The  operators  are 
assigned  to  columns  from  right  to  left  starting  with  the  equation's 
overall  output  operator  in  the  rightmost  column.  Its  input  opera- 
tors are  assigned  to  the  next  column  to  the  left  and  their  input 
operators  are  assigned  likewise.  This  process  is  continued  until  all 
operators  are  assigned  columns.  Array  OPOR  Is  used  to  keep  track  of 
this  process.  When  all  operators  are  assigned  columns,  then  Ol'K 
(1,  3)  la  defined  for  each  operator. 

OPR  (i,  A)  Is  defined  by  starting  with  the  first,  or  top,  opera- 
tor In  each  column  (right  to  left)  ami  locating  the  operator  if 
possible  so  that  the  connection  from  this  operator  to  the  one  in  the 
next  column  is  straight.  This  Is  not  always  possible  due  to  gate 
overlaps  so  the  operator  is  moved  down  the  column  until  no  overlap 
occurs.  As  each  column  is  finished,  the  next  column  to  the  left  is 
done.  This  is  continued  until  all  columns  are  finished. 

The  final  step  In  BOLD  Is  to  load  the  printout  array  PAGE  with 
the  operators,  labels,  and  connection  lines.  This  is  a straight- 
forward process  once  the  location  of  each  operator  is  defined.  In 
order  to  reduce  computer  memory  requirements,  10  characters  per  word 
are  packed  into  PAGE.  Present  dimensions  of  PAGE  are  25,200.  Op  to 
eight  pages  can  be  used  for  diagram  output  with  up  to  200  lines  per 
diagram  and  up  to  250  characters  per  line. 


BOLD  DIAGRAM 


Figure  3 


11 


CONCLUSION 


Program  BOLD  has  successfully  demonstrated  the  feasibility 
of  producing  automated  logic  diagrams  which  are  of  equal  quality 
to  manual  diagrams,  which  are  less  expensive  than  manual  diagrams, 
and  which  can  be  obtained  in  much  less  time.  Further  development 
is  required  prior  to  production  use  of  BOLD: 

- to  develop  logic  for  utilization  of  plotter  drawn  diagrams 

- to  better  define  array  dimensions  to  enable  diagramming  of 
all  equations  at  minimum  program  memory  size 

- to  add  the  capability  to  draw  feedback  lines 
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®R0GtvAM  SOLO  (INPUT  = /80  , OUT  PUT =400B , T APE4=4  0 0B, T APE7*40 OB , 
A TAPE5=INPUT,T APE6=0UTPUT) 


THIS  PROGRAM  (8-1  (ONE)  LOGIC  DIAGRAMMED  DRAWS  LOGIC  DIAGRAMS,  AS 
« • • « 

selected  by  input,  from  ri*s  emux  master  file  for  one  sioe  of  one 

AIRPLANE  EFFECTIVITY.  THIS  VERSION  OUTPUTS  LINE  PRINTER  DIAGRAMS 
ONLY. 

FILE  DESCRIPTIONS! 


TA  PEG  EMUX  MASTER  FILE  - EQUATIONS  FILE  (SEQUENTIAL) 

TAPES  INPUT  DATA  CAROS 

TAPE6  OUTPUT  FILE 

TA  PE  T EMUX  MASTER  FILE  - RANDOM  ACCESS  (SEE  PROGRAM  KANOOM 

DESCRIPTION) 

INPjr  OESCRIPTICM 

THE  ONLY  SET  OF  INFUT  DATA  CAROS  IS  IN  LIST  DIRECTED  FORMAT  AND 
IS  USFO  TO  SELECT  THE  SYSTEM,  SUBSYSTEM,  OR  SUB-SUBSYSTEM  TO 
BE  DIAGRAMMED.  THEY  ARE  SELECTED  BY  SURROUNOING  THE  FOUR  DIGIT 
OCNEE  NUMBERS  WITH  QUOTES,  BY  USING  COMMAS  AS  SEPARATORS,  AND 
8Y  PLACING  A SLASH  AT  THE  END.  SINGLE  EQUATIONS  CAN  ALSO  BE 
SELECTED  BY  INPUTTING  THE  SIGNAL  DESIGNATOR  OF  THE  DESIRED 
EQUATIONS  AS  DESCRIBED  ABOVE.  U3  TO  40  CAN  BE  SELECTED. 

AN  ADDITIONAL  FORTY-FIRST  VALUE,  IF  INPUT  AS  NON-ZERO,  WILL 
FORCE  THE  ERROR  PRINTOUT  FOR  AlL  DIAGRAMS. 

OUTOJT  DESCRIPTION! 

ONE  OR  MORE  PAGES  ARE  OUTPUT  FOR  EACH  EQUATION.  IN  AOUITIUN  TO 
HE  LOGIC  DIAGRAM,  THE  PRINTOUT  WILL  CONTAIN  THE  EQUATION  AS 
RE  AO  FROM  THE  MASTER  FIcE,  A TAB.E  OF  SIGNAL  DESCRIPTION 
INFORMATION  FOR  THE  SIGNALS  IN  THE  EQUATION,  ANO  A HEADING  WITH 
SYSTEM,  SUBSYSTEM,  AND  SUB-SUBSYSTEM  NAMES. 

VARIABLE  DESCRIPTIONS  AND  LIMITS! 

NA«E  SIZE 

CATO  30 


CM!  p 1 

CO.  12 


DESCRIPTION 

LAST  CARD  RE AO  FROM  TAPE4.  ONE  CHARACTER  PER  WORD, 
LEFT  JUSTIFIED  WITH  BLANK  FILL. 

CURRENT  CHARACTER  FRDM  EQN1 • **SEE  NOTE  10 

NEXT  RCW  INDICATOR  FDR  EACH  OF  THE  OPERATOR 
COLUMNS  IN  THE  PRINTDUT • COL(I)  CORRESPONDS  TO 
THE  ITP  COLUMN  FROM  THE  RIGHT  OF  THE  OIAGRAM. 

COL(l)  IS  THE  OUTPUT  OPERATOR  COLUMN,  COL(2)  IS  THE 
COLUHN  OF  OPERATORS  THAT  INPUT  TO  THE  OUTPUT 
OPERATOR,  ...  ,COL  (NCOL ) IS  FOR  THE  LEFTMOST 
COLUMN. 


OAT  A 400  STORAGE  ARRAY  FOR  TA’ET  RECORO  REA O BY  GETOAT. 
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fHtf  fifll  I®  BUST  QUAltTl  HhiHWlUl 
fljuN  CK>fY  furnished  toddo  — *" 


EFr 

1 

EFFECTIVITY  COOE  FOR  EQUATION.  TWO  CHARACTERS  LEFT 
JUSTIFIED  NITH  BLANK  FILL. 

sqm 

4460 

CURRENT  EQUATION  BEING  DIAGRAMMED.  ONE  CHARACTER 
PER  HORO  LEFT  JUSTIFIED  HITH  BLANK  FILL. 

••SEE  NOTE  1 

Ea*i 

960 

INTERNAL  FORM  OF  EQN.  **SEE  NOTES  2,10 

EQJ ATN 

SI  ZEE 

ARRAY  FROM  TAPE7  REC3RO  THAT  CONTAINS  THE  EQUATION 
CARO  COUNT  FOR  EACH  EQUATION  ON  THE  FILE. 

ER? 

6 

ERROR  REDEFINITION  ARRAY  USED  BY  SYSTEMC. 

♦♦SEE  NOTE  9 

IF* 

1 

LEFT  END  OF  CONNECTION  COLUMN. 

IT 

1 

EQNKII)  IS  LEFTMOST  £ NO  OF  LOGIC  TO  BE  PARSED. 

IK 

1 

X OF  (X, Y)  COORDINATES  OF  VERTICAL  SECTION  OF 
CONNECTION  PATH. 

IN)  EXM 

SIZEH 

MASTER  INDEX  ARRAY  FOR  TAPE7. 

IN’  UTS 

3f 

I NOICES  OF  VARIABLES  AND  OPERATORS  FROM  EQNl 

BEFORE  ASSIGNMENT  TO  OPR.  ••SEE  NOTES  8,10 

isrs 

1 

INPUT  SIGNAL  TO  FORCE  ERROR  PRINTOUT. 

IT  3 

1 

RIGHT  END  OF  CONNECTION  COLUHN. 

JFM 

1 

V OF  (X, Yt  COORDINATES  OF  LEFT  QASH  OF  OPERATOR 
OUTPUT. 

JJ 

1 

EQNl(JJ)  IS  RIGHTMOST  ENO  OF  LOGIC  TO  BE  PARSED. 

JT3 

1 

Y OF  (X,  Y)  COORDINATES  OF  RIGHT  DASH  OF  OPERATOR 
INPUT. 

KEf 

1 

INOEX  KEY  (NAME  TYPE)  USED  TO  GET  RECOROS  FROM 
TAPE7. 

LO°R 

1 

LAST  OPERATOR  TYPE  ENCOUNTERED.  ^SEE  NOTE  10 

NO)  L 

1 

NUMBER  OF  COLUMNS  IN  THE  DIAGRAM  (ENTRIES  IN  COL). 

NET  N 

1 

NUMBER  OF  CHARACTERS  IN  EQN. 

NE3N1 

1 

NUMBER  OF  ENTRIES  IN  EQNl. 

NIN 

1 

NUMBER  OF  ENTRIES  IN  INPUTS. 

N3’R 

1 

NUMBER  OF  OPERATORS  IN  OPR. 
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Nl»iR  1 NUMBER  OF  PARENTHESIS  SETS  IN  PAR. 

Ny MR  1 NUMBER  OF  NAMES  IN  VSR. 

OPDR  120,3  OPERATOR  ORANING  ORDER  INFORMATION.  ORDER  IN  OPOR 
DETERMINES  ORANING  ORDER  FROM  RIGHT  TO  LEFT  AND  TOP 
TO  BOTTOM. 

OPOR (1,1)  OPR  INOEX  OF  OPERATOR  IN  THIS 
POSITION. 

1 OPOR (1,21  OPERATOR  LEVEL  (*1  IS  OUTPUT 

OPERATOR,  *2  IS  INPUT  OPERATOR  TO 
LEVEL  1,  *3  IS  INPUT  OPERATORS  TO 
LEVEL  2,  ...). 

OPOR (1,3)  INOEX  IN  OPR  OF  OPERATOR  WHICH 
THIS  OPERATOR  INPUTS  TO. 

••SEE  NOTE  7 

OO?  120,34  ARRAY  OF  OPERATOR  INFORMATION.  FOR  OPERATOR  I, 

VALUES  IN  OPR  ARE  DESCRIBED  BELOW. 

OPR( 1,1)  OPERATOR  TYPE  (2003,  2004,  2005, 

2010,  2011). 

OPR( 1,2)  NUMBER  OF  INPUTS  TO  THIS 
OPERATOR. 

OPR (1,3)  IS  X OF  (X, Y)  COORDINATES  OF 
UPSER  LEFT  CORNER  OF  OPERATOR 
BOX. 

OPR (1,4)  IS  Y OF  (X,V)  COORDINATES  OF 
UP»ER  LEFT  CORNER  OF  OPERATOR 
BOX. 

OPR (1,5)  OPERATOR  INPUT  INDICES.  NEGATIVE 

THROUGH  VALUES  ARE  PRIMEO  INPUTS. 

OPR( 1,34)  ABSOLUTE  VALUES  BETWEEN  1 AND 

1090  ARE  VARIABLE  INOICES  IN  VAR. 
ABSOLUTE  VALUES  OVER  1000  ARE 
1000  PLUS  OPERATOR  INOICES  IN 
OPR. 

••SEE  NOTE  5 

PAS E 25,200  ARRAY  CONTAINING  THE  DIAGRAM  TO  BE  PRINTEO. 

••SEE  NOTE  6 

PA*  50,3  PARENTHESES  NESTING  INFORMATION.  FOR  PARENTHESIS 

SET  I,  THE  FOLLOWING  APPLIES. 

1.1  IS  THE  NESTING  LEVEL  (1  IS  INNERMOST). 

1.2  IS  THE  SUBSCRIPT  OF  THE  LEFT 
PARENTHESIS  IN  EQN1. 

1.3  IS  THE  SUBSCRIPT  OF  THE  RIGHT 
PARENTHESIS  IN  EQN1. 

••SEE  NOTE  4 

SI3F  1 ENUX  SECTION  CODE  (LEFT  OR  RIGHT). 

SI!  EE  1 MAX  NUMBER  OF  EQUATIONS  POSSIBLE  - SET  BY  GETOAT. 

SI7EM  1 DIMENSION  OF  INDEXM. 
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THIS  M89  W last  QUALITY  PRACttOBIX 

FROM  COPY  furnished  TO  DDC 

SI?  ER  1 M0R3  COUNT  OF  LAST  RECORO  READ  9Y  GETOAT. 

SI?  M 1 MAXIMUM  NUMBER  OF  KEOORDS  POSSIBLE  ON  TAPE7 • 

SY?  4C  INPUT  ARRAY  CONTAINING  OCNEE  4 DIGIT  NUMBERS  FOR 

THE  SYSTEMS,  SUBSYSTEMS,  OR  SUB-SUBSYSTEMS  TO  BE 
DIAGRAMMED. 


VAR  100,2  V AR  ( I , 1 ) IS  THE  ARRAY  OF  VARIABLE  NAMES  OR  TIME 


DELAY  DURATIONS  IN  THE  PRESENT  EQUATION.  ONE  EIGHT 
CHARACTER  NAME  PER  WORD,  LEFT  JUSTIFIED  WITH  BLANK 
FILL.  VAR(1 ,2)  INDICATES  (WITH  TWO  LETTER  LABELS) 
WHICH  VARQ,1>  NAMES  WERE  REPLACED  WITH  TWO  LETTER 
LABELS  IN  THE  DIAGRAM.  **SEE  NOTE  3 

A,R,C,  1 MISCELLANEOUS  TEMPORARIES  SET  BY  BOlO. 

0,T ,J, 

K»-  ,M, 

N,R  .S, 

T 


E , r , 6,  1 MISCELLANEOUS  TEMPORARIES  SET  BY  PARSE. 

H,  J , V,  1 MISCELLANEOUS  TEMPORARIES  SET  BY  READ4. 

W,< 

Y,?,KK,  1 MISCELLANEOUS  TEMPORARIES  SET  BY  PACK. 

LL,  MM, 

NN 

A CHANGE  IN  DIMENSION  OF  ANY  ARRAY  IN  .ABELED  COMMON  WOULO  AFFECT 
AS  A MINIMUM  - 

- ALL  LABELEO  COMMON  STATEMENTS 


- VARIABLE  DESCRIPTIONS  IN  BOLD 

- DATA  STATEMENTS  IN  BOLD 


REINITIALIZATION  STATEMENTS  IN  BOLD 


H RITE  ANO  FORMAT  STATEMENTS  IN  ECHO  ANO  ERROR 


THE  -OLL OWING  TABLE  SHOWS  THOSE  ADDITIONAL  SUBROUTINES  WHICH  MAY  BE 
AFFEDTEP  BY  A DIMENSION  CHANGE  - 


ARRAf  -—ROUTINES 

BOLO  PACK  PARSE  READ4 
CART  X 

COL  X 

EQN  X X 

EON1  X X 

ERR  X 

INPUTS  X 

OPOR  X 

OPR  X X 

PAGE  X X 

PAR  X 


oonooooooooooonnooooooooooooooonooooooonooooonooooooooo 


FOR  f HE  FOLLOWING  NOTES  - 

I IS  THE  MAXIMUM  NUMBER  OF  PARENTHESIS  PAIRS  PER  EQUATION 
J IS  THE  HAXIMUH  NUMBER  OF  UNIQUE  OPERATORS  PER  EQUATION 
K IS  THE  HAXIMUH  NUMBER  OF  INPUTS  PER  DPERATOR 
L IS  THE  MAXIMUM  NUMBER  OF  LINES  PER  DIAGRAM 
H IS  THE  MAXIMUM  NUMBER  OF  UNIQUE  VARIABLES  PER  EQUATION 
THE  PRESENT  VALUES  ARE  - 1=50,  J*120,  < = 30,  L»20Q,  M*100 

THE  FOLLOWING  TABLE  SHONS  WHICH  ARRAYS  ARE  A FUNCTION  OF  I,U,K,L,M 

ARRAY  I J K L M 


EON  XX  X 

EQN1  XX  X 

INPUTS  X 

OPOR  X 

OPR  X X 

PAGE  X 

«>AR  X 

VAR  X 


NOTE  OESCRIPT ION 

l CINENSI ON  DETERMINED  BY  2I*3J*40M. 

I CIMENSION  UCTERMINEO  BY  2I*3J*5M. 

J DIMENSION  DETERMINED  BY  M. 

. DIMENSION  0ETERM1N60  BY  1,3. 

5 DIMENSION  DETERMINED  BY  J,K*L. 

t DIMENSION  DETERMINED  BY  25, L. 

7 DIMENSION  DETERMINED  BY  J,3. 

8 DIMENSION  DETERMINED  BY  K. 

o OIMENSICN  DETERMINED  BY  COC  SYSTEM  ERROR  RESET  UiAGE. 

11  AN  INTERNAL  INTEGER  CODE  IS  USEO  TO  STORE  ANO  ANALYZE 


thf  equations. 

ALL  OR  PART 

OF  THIS  CODE  IS  USED  BY 

OPR,  EQNl,  AND 

CHAR.  THE  CODE  IS  SHDWN  8EL0W. 

VALUE 

FROM  TO 

USEO  BY 

DESCRIPTION 

-20  06 

EQNl 

COOE  FOR  NOT  EQUAL. 

-20  00  -1 

OPR, INPUTS 

PRIMED  OPERATOR  INPUT. 
ABSOLUTE  VALUE  IS 

1000  ♦ INOEX  OF 

OPERATOR  IN  OPR. 

-10  00  -1 

OPR, INPUTS 

PRIMED  VARIABLE  INOEX. 
ABSOLUTE  VALUE  IS  INOEX 
OF  VARIABLE  IN  VAR. 

1 1000 

OPR, CHAR, 
EQNl, INPUTS 

INOEX  OF  VARIABLE  IN 
VAR. 

1001  2000 

OPR, CHAR, 
EQNl, INPUTS 

VALUE  - 1000  IS  INDEX 

OF  OPERATOR  IN  OPR. 
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2001 

2002 

EQN1 

COOE 

FOR 

( 

AND 

) 

2003 

2009 

EQN1, CHAR, OPR 

COOE 

FOR 

• 

AND 

d ANO 

2006 

EQN1 

COOE 

FOR 

S 

200  7 

EQN1, CHAR, OPR 

COOE 

FOR 

t 

2010 

2011 

EQN1, CHAR, OPR 

COOE 

FOR 

$ 

A NO 

? 

9999 

EQN1 , CHAR 

COOE 

FOR 

REPLACED 

CHARACTER  - SKIPPED 
BY  LOGIC. 


VALUES  NOT  LISTED  ARE  NOT  VALIO  COOES. 

NOTE  SHOLLO  BE  HAOE  OF  THE  FACT  THAT  INCREASING  J 
OVER  1000  MILL  AFFECT  THE  INTERNAL  INTEGER 
COOE  AS  USEO  BY  BOLO  AND  PARSE.  THIS  SHOULD  BE 
AVOIOEO. 


SUBROUTINE  DESCRIPTIONS! 

NA«E  DESCRIPTION 

E2RO  DIAGNOSTIC  PRINTOUT  ROUTINE  CALLED  BY  THE  OPERATING  SYSTEM 
WHEN  FATAL  EXECUTION  ERRORS  OCCUR. 

ER?  OR  DIAGNOSTIC  PRINTOUT  ROUTINE  CALLEO  MHEN  SYNTAX  ERRORS  ARE 

FOUND  AND  MHEN  LIMITS  MOULD  3E  EXCEEDED.  EQUATION  INVOLVED 
IS  SKIPPEO.  ARGUMENT  1 TELLS  ERROR  MH AT  THE  PROBLEM  IS.  A 
NON-STANDARD  RETURN  IS  USED  10  TERMINATE  PROCESSING  OF 
CURRENT  EQUATION. 

GET  OAT  ROUTINE  TO  GET  RANOOM  ACCESS  OATA  FROM  TAPE7 . IF  DATA 

CANNOT  BE  FOUNO  THEN  QUESTION  HARivS  ARE  PUT  IN  ARRAY  DATA. 

INI T GENERAL  INITIALIZATION  SUBROUTINE.  ARGUMENT  1 IS  THE 

STARTING  LOCATION,  ARGUMENT  2 IS  THE  NUMBER  OF  LOCATIONS 
TO  BE  INITIALIZED,  ARGUMENT  3 IS  THE  INITIALIZATION  VALUE. 

LA3L  THIS  SUBROUTINE  PROVIDES  CHARACTER  (LETTER  OR  NUMBER) 

OUTPUT  FOR  USE  IN  LABELING  OPERATOR  INPUT  ANO/OR  OUTPUT 
LINES.  ARGUMENT  1 IS  A SIGNAL  TO  OUTPUT  ALPHABETIC 
OR  NUMERIC  CHARACTERS.  ARGUMENT  2 IS  THE  VALUE  TO  BE 
USEO  TO  SELECT  THE  PROPER  OUTPUT  CHARACTERS.  ARGUMENT  3 
IS  THE  OUTPUT  TNO  CHARACTERS  - LEFT  JUSTIFIED  MITH  BLANK 
FILL.  NON-STANOARO  RETURN  USEO  AS  ABOVE. 

PACK  THIS  SUBROUTINE  PACKS  CHARACTERS  INTO  ARRAY  PAGE 

(10  CHARACTERS  PER  MORO)  For  STORAGE  OF  THE  DIAGRAM.  AN 
ENTRY  POINT  (UNPACK)  PROVIOES  THE  CAPABILITY  TO  GET 
CHARACTERS  OUT  OF  PAGE.  ARGJNENT  1 IS  THE  CHARACTER(S)  TO 
8E  STORES  (PACK  CALL)  OR  CHARACTER(S)  TO  BE  RETRIEVED 
(UNPACK  CALL).  ARGUMENT  2 IS  THE  NUMBER  OF  CHARACTERS  IN 
ARGUMENT  1.  ARGUMENT  3 IS  Ti£  PRINTOUT  COuUMN  FOR  THE 
FIRST  CHARACTER.  ARGUMENT  4 IS  THE  PRINTOUT  RON. 

PARSE  THIS  SUBROUTINE  MILL  PERFORM  THE  ACTUAL  SEMANTIC  PARSING  OF 
EONl(II)  THROUGH  EQNl(JJ).  NON-STANOARO  RETURN  USED  MHEN 
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BUIS  fUftp  jg  , 

ERROR  OCCURS.  -W»R  COPY  ^isJ^^^CnCABtl 

RE1C4  SUBROUTINE  TO  SEARCH  TAPE4  FOR  THE  NEXT  EQUATION  TO  BE 

OIAGRAYMEO.  NON-STANQARO  RETURN  UiEO  WHEN  ERROR  OCCURS  OR 
FOR  NORMAL  JOe  TERMINATION. 


NRITTN  BY  R.B.  BERGER  ASO/YHEJ  9 SE3T  75  FOR  COC  6600  FORTRAN 
EXTE'JOEO  VERSION  4.4. 


REVISED  1 AUGUST  1977 
COMMON  /INFO/ 

A CA  RO(BO) ,ERR< 6)  , SYS (4  0) , ISYS,EQN(4460) , PAGE ( 25 , 200 1 , VAR < 10 0 , 2) , 
B EFF  , SlOEf 

C COL<12),ECINl(9bO>  , INPUTS  (30)  ,OPOR ( 120 , 3)  , OPR <120, 34)  ,PAR<50,3>, 
O CHAR,II,IFM,IK,  ITO,JFH,JJ,JTO,LINE,LOPR,NCOL,NEQN, NEQN1, 

E NIN, NOPR, NPAR,NVAR,A , B,C , O, E ,F , G , H , I , J,  K , KK , L, LL »M , MM, N , 

F NN,P,Q,R,S,T ,U, V,N,X,Y,Z 

C) MMON  INDEX M( 16 00 11 » EQUATN (2300) ,OATA(400> »SIZEM,SIZM,SIZEE» 

A SIZER,  KEY 
IVTESER 

A Cl  RO, ERR,  SYS,  EQN,  PAGE,  VAR, EFF,  SHE,  COL,  EQNl,OPOR,  OPR, PAR,  CHAR, 

B A,  B,C,0,E  ,F,G,H,P,Q,R,S,T,U,V,H,X,Y,Z 
II T EGER  EQUATN,OAT A, SIZEM,SIZM, SIZE-, SIZER 
EXTERNAL  ECHO 

OITA  CARO, SYS, EQN, PAGE,  VAR, EFF, SIDE  / 9782MH  / 

OAT*  ERR  / 0,0, 0,-1, -1,-1  / 

OITA  1SYS, COL, EQN1, INPUTS, OPOR, OPR, PAR, CHAR, II, IFM, IK, ITO,JFN,JU, 

A JT  0,LINE,LCPR,NCOL,NEQN, NEQN1 , NIN, NOPR, NPAR, NVAR, A ,B,C,D,E, 

B F, G,H, I, J, K, K K, L, LL, M, MH, N, NN, P, Q, R, S, T, U, V, H, X,Y,Z 

C / 5639*0  / 

SETUP  SYSTEM  ERROR  RECOVERY  LOGIC 
Cl LL  RECOVR(ECHO,77B,0) 

RE  DEFINE  ERROR  10  4 (CANNOT  FIND  RECORO  IN  TAPE7)  AS  NON-FATAL 
Cl  LL  SYSTEMC  ( 104  ,ERR) 

OPEN  T APE 7 ANO  SET  MASTER  INDEX 
Cl  LL  GETOAT ( 4H0PEN) 

POSITION  TAPE 4 INITIALLY 
RF  NIN  D 4 

RE  AO  SYSTEMS,  SUBSYSTEMS,  OR  SUB-SUBSYSTEMS  TO  BE  OIAGRAMMEO 

RE  AO<  5 , *1  SYS , IS  YS 
Ic  (SYS(l)  .NE.1H  ) GO  TO  20 
WRITE  (6,3200) 
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RE  SET  PART  OF  COMMON  INFO  BETWEEN  EACH  EQUATION 

20  C'LL  INIT (EQN, 9662, 1H  ) 

C\ LL  1NIT (COL, 12 ,- 9999) 

CVLL  INIT  (EQNl, 5626,  0) 

Gi  T NEXT  EQUATION  TO  BE  DIAGRA1M£0  - STORE  IN  EQN 
CVLL  RE  AC** , RETURNSI20, 999) 

Srl  EQN1  A NO  VAR 

TL  ST  FOR  OPERATOR  - IF  FCUNJ,  THEN  PUT  IN  PROPER  CODE  NUHBER  - 
( IS  29  C 1 
) IS  2002 

• IS  20C3 
0 IS  200  4 

♦ IS  20  OS 

* IS  2006 

• IS  2007 
f IS  2010 
T IS  2011 

IF  EQN (I)  IS  NOT  AN  OPERATOR,  THEN  IT  IS  EITHER  A VARIABLE 
NAME  OR  A TIME  DELAY  DURATION.  VARIABLES  AND  TINE  OELAVS  MILL 
BE  REPLACED  BY  THEIR  LOCATION  IN  VAR. 

I'-  1 

30  j=o 

Ir  (E3N(I)  .EQ.  1M< ) J=  2 00 1 
I- (EQN(I) .EQ.1H) ) J= 2002 
tr  <E2MI»  .EQ.1HM  J=  200  3 
Ir (EQN(I) .EQ. IMS ) J*20Q4 
I r ( E3  N ( I ) . EQ  . 1H ♦ I J*  200  5 
Ir  (E3NCI) .EQ.1H*)  Js  2006 
I r (EQN ( I ) .EQ • 1H ' ) J- 2007 
1-  <EON<I» .EQ.1HSI  J=2010 
Ir  (EQN(I)  .EQ.1H?)  J=201 1 
Ir  ( J.  EQ.0)  GO  TO  40 

OPERATOR  FOUND  - SET  EQN  1 (NEQN1)  AND  INCREMENT  I 
N*  QN1  = NEQNl* 1 

IF  (NEQN1.GT.960)  CALL  ERROR(2>,  RETJRNS(20) 

E3NKNEQN1)*  J 
1*1 

Ie  < I- NEQN ) 30,30,140 

EITHER  VARIABLE  NAHE  OR  TIME  DELAY  DURATION  STARTING  AT  EQN(I) 
T:ST  FQN ( I *4 ) TO  FIGURE  OUT  WHICH 

40  Ir  (EQN(I«4».NE.1H-  .AND.  EQN  < I ♦ 4> .NE • 1H>  .AND. 

A <FQN(I*4> .LT.1HA  .OR.  EQN  (1*4) . GT . 1HZ> ) GO  TO  70 
C 
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VARIABLE  NAME  IN  EON(I)  THROUGH  EQN(I*7)  - PACK  INTO  VAR(NVAR) 
IF  IT  IS  A ME  H NAME  ANO  THEN  SET  EQNl(NEQNi)  ANO  INCREMENT  I 

EH  COOL  (8,  200  0,K)  EQN  ( I)  , EQN  ( IU ) , EQM  ( I *2 ) , EQN(I*J) , EQN  ( 1*4  > , 

A E3N(I*5),£QN(I»6),EQN<I*7) 

03  5C  J* 1 ,NV  AR 

50  Ir  (K.  EQ.VAR!J,1)  ) GO  TO  60 

Tc (NV  AR.EQ.l 00)  CALL  ERROR! II , RETURNS (20) 

NV  AR*  NVAR»1 
VI  R(N  VAR,  i)*K 
J*  NVAR 

60  N-ONl  = NEQNltl 

Ic  ( It  QN1 . GT.  960 ) CALL  ERROR  (2)  , RETJRNS(20) 

E3  Nl(  NEON  1)  * J 
I*  !♦§ 

Ir  ( I- NEQN)  30*30,140 

TIHE  DELAY  OURATICN  IN  EQN(I)  THROUGH  EQN(I*3>  - PACK  INTO 
VAR (NVAR)  IF  IT  IS  A NEM  VALUE  AMO  THEN  SET  EQNl(NEQNl)  AND 
INCREMENT  I 

70  EVC03E(4,2000*KI  EQN  ( I)  * EQN  ( 1 + 1)  * EQN  (I+2),EQN(I»3) 

0)  8P  J*1 * NV  AR 

80  I" (K. EQ. VARC J,l) ) GO  TO  90 

Ie  ( NV  AR.  EQ.  1 0 0)  CALL  ERROR!  1)  * RETURNS(2J) 

NV  AR=  NVAR*1 
V'  R (N  VAR, 1) *K 
J=  NVA R 

90  m-:qni=neqni*  i 

Ir  (NEQN1.GT.  960)  CALL  ERROR  (2)  * RETJRNS(20) 

E5NKNEQN1I*  J 
Is  I *4 

IF ( I- NEQN)  30,30,140 

SCAN  EQN1  FOR  PARENTHESES  TO  DETERMINE  PROPER  PARSING  OF 
EQUATION  - SET  PAR  ANO  NPAR 


140  J : K*l.  * 0 

01  17  C I = 1,N£CN1 

LI  OK  FOR  LEFT  OR  RIGHT  PARENTHESIS 

l~  (EQNKI1.EQ.2002*  GO  TO  150 
IF  (EQNl(I)  .NE.2001)  GO  TO  170 

LEFT  PARENTHESIS  FCUNO  - SAVE  SUBSCRIPT  OF  MOST  CURRENT  LEFT 
PARENTHESIS  IN  J,  INCREMENT  K (THE  NESTING  LEVEL),  SET  L (THE 
MAXIMUM  NESTING  LEVEL!,  AND  SET  »AR(J,1)  AND  PAR(J,2) 

l-  (NEAR. FQ. 50)  CALL  E FROR! 3 ) , RETURMS(20) 

N>  AR=  NPAR  *t 
Ji  NPAR 
K 1 K*1 

L : M AX t (L,  K) 

°1R(J,1)=K 
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PAR  (J  » 2)  = I 
GD  TO  170 

RIGHT  PARENTHESIS  FOUND  - SET  PAR(J,3),  OECREHENT  K (NESTING 
LEVEL)  ANO  RESET  J (RIGHTMOST  LE“T  PARENTHESIS  FOUNO  WITHOUT 
A RIGHT  PARENTHESIS) 

ISO  Ir(J.  EQ.O)  CALL  ERROR  (<»)  , RETURNS  ( 23  ) 

<=  K-l 

P»  R(J,3)  = I 
160  J-J-l 

IF  (J.  EQ.O  GO  TO  170 
Ir  <PAR<  j,3).NE.O)  GO  TO  160 
170  CD  NTI  NUE 

CH  ECK  FOR  UNHATCHEC  PARENTHESES 

Ir(J.NE.0)  CALL  ERROR  (8 ) , KETURNS(20) 

Ic(W»AR-i>  2LO,220,180 

RESET  NESTING  LEVEL  AS  STOREO  IN  PAR  SO  THAT  INNERMOST  SET 
H»S  ONE  FOR  ITS  LEVEL 

180  CD  190  1=1, NEAR 
190  °AF  (I  ,1)=L-PAR(I,1)*1 

ARRANGE  PAR  FROM  LCW  TO  HIGH  LtVEL  ANO,  FOR  EACH  LEVEL,  FROM 
LEFT  TO  RIGHT  POSITION  OF  LEFT  PARENTHESIS. 

K=  NPA  R-l 
CD  21  i 1=1, K 
L=  m 

CD  210  J=L,NFAR 

I-  (PAR(I,l)*  1000  0*PAR(I,2).LE.PAR(J,l)*10000fPAR(J,2)>  GO  TO  210 
CD  20  0 M=  1,3 
N = PAR (I,M) 

®AR  a ,M)=PAR(J,M) 

200  °5R(J,M)=N 
210  CD  NTINUE 

FOR  EACH  PARENTHESIS  SET,  CALL  PARSE  TO  ANALYZE  LOGIC  IN 
EQN1 ( II)  THROUGH  EQNKJJ)  ANO  THEN  WIPE  OUT  PARENTHESES. 

220  CD  23  0 I = 1,NPAR 
It  = PA  F ( 1 , 2)  ♦ 1 
JJ=PAR(I,3)-1 
CALL  PARSE,  RETURNS  (20) 

230  E3Ni(II-l)*EQNi( JJ»1>=9999 

NO  MJFt  PARENTHESES  - PARSE  ENTIRE  EQUATION 

2 <»0  ir=3 

JJ  =NEONl 

C»  LL  PARSE,  RETURNS  ( ) 
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THE  EQUATION  HAS  BEEN  PARSED  - NOW  NEED  TO  DETERMINE  WHERE  THE 
ITERATORS  WILL  BE  LOCATED  ON  THE  PAGE  PRINTOUT  TO  MAKE  THE  MOST 
READABLE  OIAGRAM  (MINIMUM  CROSSOVER  OF  OPERATOR  CONNECTIONS, 
MINIMUM  NUMBER  OF  UNCONNECTED  OPERATORS,  MAXIMUM  USE  OF 
STRAIGHT  LINE  CONNECTIONS,  ALL  CONNECTIONS  GOING  LEFT  TO  RIGHT) 

CHECK  FOR  ONLY  ONE  OPERATOR  IN  THE  EQUATION 

Ir  ( NOPR-1 ) 20,250,260 
250  O3  R (1 , 3)  = 10 

03R(1 ,4)=0PR  (1,2)*MOO(OPR(1 ,2)  ,2)-l 

n: 01=1 

COL  (1  )*OPRCl  ,4)*0PR<1,2)  -MOD  (OPR  ( 1,  2 ) ,2)  *3 
IE  (C0L(1)  .GT.203)  CALL  ERROR(ll),  RETURNS(20) 

SO  TO  400 

MORE  THAN  ONE  OPERATOR  IN  THE  EQJATION.  BEFORE  OPR(I,3>  AND 
OPR  ( 1, 4)  CAN  BE  SET,  NEEO  TO  DEFINE  OPOR(I,l),  0P0R(I,2),  AND 
OP  OR  (1 , 3)  • ORDER  CF  ENTRIES  IN  OPOR  IS  OIAGRAM  ORDER  STARTING 
WITH  OUTPUT  OPERATOR  ON  THE  RIGHT  ANO  WORKING  TO  THE  LEFT  - 
TOP  TO  BOTTOM  IN  EACH  COLUMN. 

OPOR  (1,1)  IS  THE  INOEX  OF  THE  OPERATOR  IN  OPR 

OPOR (1,2)  IS  THE  OPERATOR  LEVEL  (=1  FOR  OUTPUT  OPERATOR, 

=2  FOR  INPUT  OPERATORS  TO  LEVEL  1,  =3  FOR  INPUT 
OPERATORS  TO  LEVE.  2,  ...) 

OPOR  (1,3)  IS  THE  INOEX  IN  03R  OF  THE  OPERATOR  WHICH  THIS 
OPERATOR  INPUTS  TD 

260  M=  OPOR(l,  2)  = 1 
Q3  OR( 1,1) =NO  PR 

TAKE  12  PASSES  THROUGH  OPR  - ONE  FOR  EACH  LEVEL.  IF  ALL  OF 
OPOR  IS  NOT  SET,  THEN  THE  EQUATION  REQUIRES  MORE  THAN  12  COLUMN 

DO  272  L=l,12 

FIND  LEVEL  L IN  OFCR  AND  THEN  CO  TO  OPERATOR  J IN  OPR  AND 
ASSIGN  ITS  INPUT  OFEKATORS,  IF  ANY,  TO  LEVEL  L*1 

IT  2*1  J= 1 ,N  OFR 

Ir  (OP  OR  (J, 2)  .NE.L)  GO  TO  271 

S*  OPOR  <J,1> 

N-OPR  (S, 2) 

DO  27  G K=1,N 
R=  IA3S<OPR(S,K»4>) 

IF (R. LT.1000)  GO  TO  270 
NsP*l 

OPOR(  M, 1) =R- 1000 
O3 OR(  H, 2)  =L*  1 
O3  OR(  M,  3)  *S 

Ic (M. EQ.NOPR)  GO  TO  290 

270  CONTINUE 

271  CONTINUE 

272  CONTI  NUE 
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FWLEC  TO  ASSIGN  ALL  OPERATORS 
O'  LL  LRROR(fe),  RET  OR  NS ( 2 0 ) 

SFT  NCOL  SO  THAT  IT  IS  THE  MAX  LEVEL  (NUMBER  OF  COLUMNS) 

290  N: 0L=0P0R (1.2) 

T)  3C  C 1=  2, NOPR 
300  9!.  OL=hAXO  (NCOL.OPOR  (1 ,2)  * 


SFT  OPR  (I  .3)  SO  THAT  THE  OPERATORS  HILL  START  IN  COLUMNS 

lr,  3C  . SO,  TO,  90  , 110,  130,  150,  170,  190,  210  , 230  IN  ARRAY 

PI  GE 


93  31  t I = 1 , N O PR 

31 P 03F  (OPOR ( 1,1)  ,3)  =20*  (NCOL-OPOMI  ,2)»1) -10 

SM  OPR(I,A>  BY  GOING  THROUGH  OPOR  IN  OROER.  PROCEDURE  IS  TO 
PLACE  THE  OPERATOR  VERTICALLY  SO  THAT  THE  CONNECTION  PATH  IS 
STRAIGHT  - IF  POSSIBLE.  OPERATORS  ARE  PLACEO  TOP  TO  BOTTOM 
WITHIN  A GIVEN  LEVFL  - COL  IS  US-IO  TO  KEEP  TRACK  OF  THE  NEXT 
AVAILABLE  POSITION  IN  EACH  LEVEL  SO  THAT  THE  OPERATORS  WILL 
HAVE  AT  LEAST  TWO  LINES  BETWEEN  THEM. 

SlT  OPMNOPR.A)  ANC  THEN  DO  THE  REST 

O1  R (NOPR,  A)  = OFR(NOPR,2)  ♦MOO  (OFR(NOPR  ,2)  ,2)  -1 
Old  ) =OPR(NOPR,  A)  fOPR  (NOPR  ,2)  -MOD  (OPR  (N0PR,2) , 2)  *3 

SET  OPR  (I  , A)  - NEEO  TO  KNOW  ITS  INPUT  POSITION.  THE  FIRST 
EXPRESSION  IN  THE  MA X 0 IS  THE  Hi 3HEST  AVAILABLE  POSITION 
WITHOUT  OVERLAP  ANC  THE  SECOND  EXPRESSION  IS  THE  DESIREO 
VERTICAL  POSITION  FOR  A STRAIGHT  CONNECTION  PATH 

•n  33  0 1=2, NOPR 
S-  OPOR  (1,3) 

N;  OPR (S, 2) 

00  3?t  J=  1 ,N 
R-  IAfcS (OPR (S , J*A ) ) 

Ir  (R.  NE.OPOR  (I,l)*1000)  GO  TO  320 

R = R-l COO 

<=OF3MI,2) 

H=  OPR (R, 2) 

O'  R (F  , A ) = OPR  (S»4)-N-NC0(N,2)*2*J-1 

IF (CoL(K) .NE • -9999)  OPR («, A ) =MAX0 (OaR ( R, A) , COL ( K) ♦M*MOD(M, 2) -2) 
CO  L (<  ) *OPR  (R , A)  ♦M-MOO  (M , 2)  ♦ 3 
GO  TO  330 
320  CONTINUE 
330  CONTINUE 

NOW  THAT  OPR( I, A)  HAS  BEEN  SET  (BASED  ON  THE  OUTPUT 
OPERATOR  BEING  IN  LINE  1),  NEEO  TO  INSURE  THAT  THE  TOP  OF 
THE  HIGHEST  OPERATOR  IS  IN  LINE  1 

Lc T L BE  ThE  TOP  OF  THE  HIGHEST  OPERATOR 
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L * OPR (1,4)-0PR(1,2)-H00( OPR (1,2) ,2) *2 
0 ) 340  I>2,N0PR 

34  0 L-  MIMO(L,OPR(I#  4 >-0PR (I , 2) -MOO (0PR(I ,2)  ,2)  *2) 

NOM  NEEO  TO  RESET  0PR(I,4)  AND  CDL  < I> 

03  35 0 I*1,NC0L 
350  Old)  ■COL  (I  ) -L*  1 
03  3b  0 I>1,N0PR 
360  O’R (I ,4)*0PR(I,4)-L»1 

CHECK  COL  TO  SEE  THAT  THE  0IAGRAM  HILL  FIT  IN  ARRAY  PAGE 
03  39b  Is  1,NC0L 

390  I" (COL(I) .GT .203)  CALL  ERR0R(11)»  RETURNS(20) 

LOAD  OPERATORS  INTC  PAGE  AND  THEM  CONNECT  THEM 

40  0 0 3 59  0 N«1,N0PR 
A=  OPR (N,2) 

J 3 OPR <N,  3) 

<*  OPR  (N,  4) 

L=K-4-M0D(A,2)»2 
3= K *A -HOD (A, 2 ) 

OUTPUT  THE  SOX  (EVERYTHING  EXCEPT  INPUT  OASHES,  PRIMES,  AND 
LEFTMOST  ASTERISKS) 

CALL  FACK(2H*%2,J*1,K)  , RE  TURNS (20) 

Ir <OPR(N,1).EQ.2003>  CALL  P ACK (3HA*- , 3 , J*l, K*l) , RETURNS(20> 
l-  ( CPMN,  1)  .EQ.2004)  CALL  PACK ( 3HX*-, 3. J»l, K»l> , RETURNS  ( 20 ) 

Ir (CPR(N,1).EQ.2Q05)  CALL  PACK( 3HO* - , 3 , J*l, K*1 ) , RETURNS(20) 

I-  (CPR(N,  D.EQ.2010)  CALL  PACK  < 3H1*- , 3 , J*l,  K*l>  , RETURNS(20) 

I~ (OPK(N,l>.EQ. 20111  CALL  PACK (3H2*-» 3, J*l, K*l) » RETUKNS(20) 

CALL  PACK(2H**,2,J»1,K*2>,  RETURNS(20) 

SET  THE  INFUT  OASHES,  PRIMES,  AN3  LEFTMOST  ASTERISKS 

0 3 42  0 M=L ,D 
S=  MOP  (M-L *1,  2» 

Ic (S. EQ. 0 ) CALL  PACK ( 2H  *,2,J-1,M),  RETURNS(20) 

Ir  <S. NE.O  .ANO.  OPR(N,(M-L*2)/2>4> .GE.O)  CALL  PACKT 2H-* ,2 , J-l ,M ) , 
A RE  T URNS  (2  0) 

420  IF  (S.  NE.O  .AND.  OPR (N , ( M-l» 2>/2 *4) . - T . 0 > CALL  PACK(  2H-0 ,2  , J-l , H ) , 
A RFT  URNS (20 ) 

CONNECT  THE  OPERATORS  ANO  LABEL  THE  INPUTS  FROM  VARIABLES. 

THE  ONLY  LINE  CONNECTION  PATH  TO  BE  ATTEMPTED  WILL  BE 
OVER-UP (004N) -OVER.  IF  THERE  IS  NO  CLEAR  PATH  AVAILABLE,  THE 
OPERATOR  OUTPUT  ANC  INPUT  HILL  BE  LABELEO  HITH  A T HO 
01  GIT  NUMBER.  OPERATOR  INPUTS  FROM  VARIABLES  HILL  BE  LABELEO 
HITH  A THO  LETTER  LABEL.  A TABLE  HILL  BE  PRINTED  TO  IDENTIFY 
T( t THO  LETTER  LABELS  HITH  RI»S  SIGNAL  DESIGNATORS. 
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LOOK  AT  EACH  INPUT  TO  EACH  OPERATOR  ANO  CONNECT  OR  LABEL 
THOSE  INPUTS  FROM  CTHER  OPERATORS.  INPUTS  FROM  VARIABLES 
HILL  BE  OONE  NEXT. 


03  59  C 1*1, A 
L=  IA3S(OPR(N,I*4)> 

IF (L. Lf  00)  GO  TO  590 
L*L-1  00b 


INPUT  IS  FROM  ANOTHER  OPERATOR  - TRY  TO  CONNECT  THE  OPERATORS 


iro*j-5 

jr  0=K-A-NOO(  A, 2)  *2*1 
IF M*0PR(L,3) *7 
Jr  M*OPR (L ,4) *1 
S = HIN0  ( JFH,  JTO) 

T = PAX  C (JFM,  JTO) 

FIND  CONNECTION  PATH 

03  530  IK  = IFM,IT 0 

SEE  IF  PATH  IS  CLEAR  FROM  ( IFM, JFM)  TO  (IK, JFM)  - BLANK  OR 
I OK 

03  49  0 M=  IFM , IK 

CALL  UNPACK! C,1,M,JFH) , RETURNS (20) 

490  Ir (C.  NE. 1H  .ANO.  C.NE.1HI)  GO  TO  590 

SEE  IF  PATH  IS  CLEAR  BETWEEN  (IK,  JFM)  ANO  (IK,  JTO)  - BLANK  OR 
08  SH  OK 


35  5T0  M=S,T 

CALL  UNPACK!  C ,1 , IK  ,M ) , RETURNS(20) 

500  IF(C.  NE.1H  .AND.  C.NE.1H-)  GO  TO  590 

SEE  IF  PATH  IS  CLEAR  BETWEEN  (IK, JTO)  ANO  (ITO,JTO)  - BLANK  OK 

03  51 C N=IK,  ITO 

C»LL  UNPACK  ( C,1,M,  JTO),  RETURNS(2u) 

510  IF  (C.  NE.  1H  ) GO  TO  58C 


CONNECTION  PATH  IS  CLEAR  BUT  BEFORE  DRAWING  NEEO  TO 
MIKE  SURE  THAT  IK  IS  IN  THE  PROPER  COLUMNS  (17-25,  37-45, 
57-65,  77-85,  97-105,  117-125,  137-145,  157-165,  177-185, 
19  7-205,  217-  225 J 


IF  ( ([ 

K, 

»GT 

. 25 

.AND. 

IK, 

.LT. 

37) 

.OR. 

(I 

K. 

> GT . 

45 

.ANO. 

IK 

.LT. 

57) 

.OR 

<r 

K. 

,GT 

. 65 

.ANO. 

IK. 

,LT. 

77) 

. OR  • 

(I 

K, 

, GT , 

85 

.AND. 

IK 

.LT. 

97) 

.OR 

a 

K, 

>GT 

.105 

.ANO. 

IK. 

,LT. 

117) 

.OR. 

(I 

K « 

, GT , 

125 

.ANO. 

IK 

.LT. 

137)  , 

.OR 

(i 

K, 

i GT 

.14  5 

.ANO. 

IK, 

>LT. 

157) 

.OR. 

(I 

K. 

> GT  a 

165 

.ANO. 

IK 

.LT. 

177) 

.OR 

a 

K, 

,GT 

.185 

.AND. 

IK, 

,LT. 

197) 

.OR. 

(I 

K. 

> GT , 

205 

.AND. 

IK 

.LT. 

217)  1 

> 

E GO  TO  580 


PATH  FOUNO  WITH  CORNERS  AT  (IK,J~M)  ANO  (IK, JTO)  - LOAO  IT  INTO 
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non  o o o o noon  o o o o o o 


PA  GE 


- USE  PERIODS  FOR  CORNER  CHARACTERS 


L* IFM-3 
0* IT3  *3 

I?  <JFH.NE.JTO)  GO  TO  530 
03  5? 0 M*L,D 

OLE  UNPACK<C,1,M,  JFM),  RETURNS<20) 

520  IFtC.EQ.lH  ) CALL  PACK<  1H-,  1,M,  JFH) , RETURNS(20) 
GO  TO  59P 
530  S*S*1 
T*  T -1 

0)  5?C  N«S,T 

550  C*  LL  PACK <1HI,1,I< ,M>  , RETURNS<20> 

03  56C  H*L»I  K 

OIL  UNPACK<C,t,M,  JFP),  RETURNS(20) 

560  I?  (C.EQ.iN  I CALL  PACK < iH-, 1,H, JFM)  , HETURNS(ZO) 
0 3 57  u H*IK,0 

OIL  UNPACMC,1,N,  JT0)  , RE  T URNS  ( 20  ) 

570  I r <C.  EQ. 1 H I CALL  PACK < IN- , 1 ,M, JTO) , RETURNS<20) 
CALL  PACKC1H.,1,  IK, JFM  , RETURNS(20) 

CALL  PACM1M. ,1, IK, JTO) , RETURNS(20) 

S3  T3  590 
5S0  CONTINUE 


NO  PATH  POSSIBLE  - INSTEA0  OF  CONNECTING  OPERATORS,  LABEL  THE 
OUTPUT  AND  INPUT  WITH  A TWO  DIGIT  NUMBER. 


CALL  LABL(2,L,3),  RETURNS<20) 

CALL  PACKIB, 2, IFH-3, JFM)  , RETURNS<20) 

CALL  PACK <B, 2 ,IT 0* 2 , JTO ) , RETURNS<20) 

590  CONTINUE 

LOOK  AT  EACH  INPUT  TO  EACH  OPERATOR  AND,  FOR  THOSE  INPUTS  FROM 
VARIABLES,  LABEL  WITH  V AR(L,1>  OR  WITH  A 2-LETTER  LABEL 

0 3 62 C N* i»N  OPR 
A*  OPP (N,2) 

J*  OPR  <N,3» 

<-  OPF  <N,4 ) 

0)  620  1*1,  A 

L< I ABS (OPRCN , 1*4 ) ) 

Ir  <L.  GT.10  00)  SO  TO  620 

INPUT  IS  EITHER  A VARIABLE  OR  >\  TIME  DELAY  DURATION  VALUE  - 
LABEL  WITH  VAR(L,1)  OR  A 2 LETTER  LABEL 

T * K-A-MOO(A,2)f2*I 

If  < OPR<N,l)  .NE.2010  .ANO.  OPR  (N , 1 ).  NE.  2011 ) .OR. 

A AN 0 ( VA R (L  , 1)  ,7777777777778) .NE . 555555555555B)  GO  TO  600 

INPUT  IS  TINE  DELAY  DURATION 

CALL  UNPACK <C  ,5, J-6, T ) , RETURNS<20> 

IF (C.  NE. 1H  ) GO  TO  610 

CALL  PACK  (VAR  <L»  II  »%»  J-5,T)  , RETURNS(20> 
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OTIS  PAGE  IS  BEST  QUALITY  PRACTXCA&Lfi 

•Clf  rr  * y'  ' ' OTOM  COPY  FURNISHED  10  JlDfl  " 

G?  TO  620 
C 

C INPUT  IS  A VARIABLE 

C 

600  I-U.EQ.IO)  CALL  UNPACK (C, 8 ,J-9,T) , RETURNS(20) 

Ic (J.  ME. 10)  CALL  UNPACK(C,9, J-10,I>,  RETURNS(20> 

Ir(C.NE.lH  ) GO  TO  610 

Cl  LL  PACK(VAR(L,i>  ,S,J-9,T)  , RETURNS(20> 

GD  TO  620 
C 

C INPUT  COULO  NOT  BE  FIT  IN  - USE  A 2 LETTER  LABEL 

C 

610  GILL  LABL  (1 ,L  ,VAR(L , 2 ) ) , RETUKNS(ZO) 

C»LL  PACK(VAR(L,2>  ,2,J-3,T)  , RETURNS(20) 

620  CONTINUE 
C 

C ADO  OUTPUT  VARIABLE 

C 

i-  OPR (NOPKf  3 ) 

<=  OPR (N0PR,4) »1 

IF (EJN1 (2).EQ.-2006)  CALL  PACK< 1HO, 1 , J*2,K) , RETURNS(20> 

SILL  PACK  (VAR  <1, 1)  ,B,  J4*,K)  , RETURNS  (20) 

C 

C OUTPUT  HE  ACER 

C 

Cl  LL  GETDAT( VAR (1,1)) 

0:  CODE ( 16  V , 2400 , OAT  A)  L,M,N 
S*rODE(10,1900,KET)  VAR(i,l> 

GILL  GETOAT  ( KEY) 

LINE'S 

WR  IT*  (6,2800  I VAR(1, 1)  ,EFF,  SIDE,  (DATA  (I)  ,1*1,7)  ,L,H,N, 

A (DATA  (I),  1*8,21) 

C 

C OUTPUT  EQUATION 

C 

J = 125 

Ic  (N:  QN.GT.124)  GO  TO  630 
LTNE=LINE*2 

WRITE  (6,2300)  (EQN(K)  , K = 1 » N EQN ) 

GD  TO  670 
630  J*J-i 

Ie  (EQN( J)  .NE.1H*  .AND.  EQN( J) .NE.1HD  .AND.  EQN ( J) .NE. 1H*  .AND. 

A EiN(J) .NE.1HS  .ANO.  EQN( J) .NE.1H7)  GO  TO  630 

LI  NES  LINE *2 

MR  ITS (6,2300)  (EQN (K) ,K*1, J ) 

640  I*JU 

Ir  (Jd24.LT.NEQN)  GO  TO  650 
LI  NE*  LINE*1 

WRITE  (6,2200)  (EQN(K)  ,K*I,NEQN) 

G)  TO  670 
650  J * J * 1 2 5 
660  J!  J*1 

IF  (E3N(J) .NE.1M*  .ANO.  EQN( J) .NE. 1HJ  .ANO.  EQN ( J) .NE. IN*  .AND. 

A E3N(J) .NE.1HS  .AND.  EQN( J) .NE. 1H7  .AND.  I.NE.JI  GO  TO  660 
LI  NE*LINE  *1 

j 
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uug  uuu  ooo 


HR ITE  (6, 2200 ) (EQN(K> ,K*I, J) 
GO  TO  6V0 


SSfttSJii?81  QUALm  PWCriCABU 

OOPT  FURNISHED  IX)  DDC 


REARRANGE  VAR  FROH  LON  TO  HIGH 

670  I * NVA  R-i 

03  680  J*2,I 
<’  J*1 

0 0 68  C l*K,NVAR 

V (V4R(J,1).LE.VAR(L,1>>  GO  TO  680 
N = VAR ( Jy 1 ) 

Vi  R (J  * 1 1*  VAR  CL*  1 0 
ViRt.  »1)*M 
1*  VAF ( J,2) 

Vi  R(J  * 2)  = VAR  <1,2 ) 

ViR(l»2)*M 
680  30  NTI  NUE 

OUTPUT  SIGNAL  DESCRIPTION  TABLE 


NR  XTE (6*7  000 1 

0s  NVA  R/2*HOO  (NVAR,  2) 

LI  NE=  LINE *B* 6 
00  7CC  1*1, B 
A*  C*D*J*K=R=S=T=  1H 
L*  1 OR  TINE  DELAY 
N*  10H  DURATION 
Ns  1 CH  (SEC) 

Ic  (AN  Cf  VAR  (I  ,1)  ,77  777777  777  7 B)  • EQ  . 5555555555558)  GO  TO  690 
CiLL  OETOAT ( VAR( I, 1) ) 

9ECODEU60,2AOO,OATA)  L , H,N 
3:COOE(20,2100,VAR(I,2))  J,K 
690  IF (8*I.GT.NVAR  .OR.  NVAR.EQ.l)  GO  TO  700 
A--10HTIHE  DELAY 
0*  ION  DURATION 
S* 10H (SEC) 

C*  VAN  (8H,1> 

Is  (AN  DIC,  7777777  77  77  7 B)  • EQ.  565565i5i  i.  *»  'B)  GO  TO  700 
CiLL  GETDAT(C) 

O-I  COO  E ( 16  0 , 2 LOO  , DAT  A)  A,0,S 
DEC OD L( 20, 2100, VAR (8+1, 2))  T,R 
700  HR ITE (6, 3100 ) J, K, VAR <1 , 1) , L,M, N,  T,R,C,A,D,S 


OUTPUT  DIAGRAH 


J=  NCOL-6 
<= A*- 9999 
DO  710  I*l,NCOL 
IF(I.LE.J)  A*MAXO(A,COL  (I)) 

710  I~((I.GE.U  .ANO.  J.GT.C)  .OR.  J.LE.O)  K*MAX0 ( K,COL ( I) ) 
CiLL  GETDAT( VAR( 1,1) ) 

0: CODE (16  0, 2 A CO,  DATA)  L , N,N 
EMC  ODEIIO, 1900, KEY)  VAR(1,1) 

CiLL  GETDAT (KEY) 


30 


ooo  r>  o o non  non  non  on 


, THIS  PAGE  IS  BEST  QUALITY  PR1CTICAIL1 

„ . ..  --  •'  raOH  COH  FUKM1SH8D  TODDC  

SEE  IF  HEADER’ PAGE  CAN  BE  USED 

Ir  (LINE*K.GT.55>  GO  TO  73C 
HilTE  (6,2700  ) 

Hi ITE (6,2900 ) ( (PAGE (I , J) , I *1,13) ,J*i,K) 

HEADER  PAGE  USED  - SEE  IF  ADDITIDNAL  PAGE(S)  NEEDED 

le  (A.  EQ.-9999)  GO  TO  740 
IF  (LI)OK.LE.51)  HRI TE (6  ,26 00) 

Hi ITE (6,2600)  HAR( 1,1) ,EFF, SIDE, (OAT A (I) ,1=1, 7) ,L ,H ,N , 

A (DATA  (I)  ,1=8,21) 

J = LINE-7 
3)  72  C 1=1, J 
720  Hi ITE (6,3400) 

HilTE  (6,3300)  ((PAGE C I , J) , I =14, 25) , J = 1 , A) 

GD  TD  740 

HEADER  PAGE  CANNOT  BE  USEO  - GO  TO  A N£H  PAGE 

730  I- (LINE. LE. 51)  HRI TE (6, 2500 ) 

H?IT£ (6,2600  ) HA R ( 1 , 1 ) , EFF , SIDE , (OATA(I) ,1=1,7) ,L  ,M , N , 

A (3 ATA(l), 1=6, 21) 

Hi ITE  (6,2700) 

Hi ITE  (6,2900)  ( (PA GE (X, J ) ,1 =1, 13) , J=1 , K) 

HEADER  PAGE  NOT  USEO  - SEE  IF  ADDITIONAL  PAGE(S)  NEEDED 

I-  (A.  EQ. -9999)  GO  TO  740 
Ir  (LINEAR. LE. 51)  HRITE  ( 6 ,26  00) 

Hi ITE (6,2600)  HA R( 1 , 1 ) ,£FF, SIDE, (DAIA(I) ,1  = 1,7) ,L,H,N, 

A O AT  A (I),  1=8,21) 

Hi ITE (6,2700) 

Hi ITE  (6,3300)  ( ( PAGE (1 , J ) , I =14 , 25) , J= 1 , A ) 

THIS  EQUATION  IS  CONE  - GO  TO  NEXT  ONE 

740  Ie (IS YS.NE.C)  CALL  ERROR  (7)  , RETURNS(20) 

GD  TD  20 

NO  HOkE  EQUATIONS 

999  CD  NTINUE 
1900  FDRMAT(A4*  LABEL*) 

2000  pDRHAT (8A1) 

2100  F)  RMA  T (2A1) 

2200  FDFMAT  (11X,124A1) 

2300  FDRMAT(/2X,124A1) 

2400  FDPHAT  (34X,2A10,  A5  ) 

2500  F) RHAT (//50X*Stf$$SS3$SlSS$S»$$$$5$$$Si$S$SiiSi*/ 

A 50X*SSfS  OIAGRAH  ON  THE  NEXT  PAGE  $S$$*/ 

B 50X*mSSSS*!SfSS$mS*$$St$StSSmSSS*> 

2600  FDRMflT (//50X*i*SSSf$Sff $SS$ *$$$$$ $$$*$${»$$*$$$*/ 

A 50X*im  Cl  AGRA  M CONTINUED  $$**•/ 

p »ox*msimmnmmsmmmmit*) 
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270P  F3RNAT<*0*I 

2600  nRMAT<*l*///23X*EMUX  BOOLEAN  EQUATION  DIAGRAM  FOR  SIGNAL  *A8 
A ♦ OF  AIRCRAFT  NO.  *A2,2X,A5*  SIDE  OF  EMUX*// 

B 1&X.7A10,  5X.2A10.A5/ 

C 15X.7A10/ 

0 15X.7A10) 

2500  F)RMAT(IX.IJAIO) 

SCO**  FORMAT  (/// 

A 13  X*  DIAGRAM  EMUX  SIGNAL  SOURCE  OR  DESTINATION* 

5 10X*  DIAGRAM  EMUX  SIGNAL  SOURCE  OK  DESTINATION*/ 

C 10 X* DESIGNATOR  SIGNAL  * 

0 A2X*OFSIGNATOR  SIGNAL*/ 

F 1GX*  DESIGNATOR* 

F *IX*  DESIGNATOR*/) 

3100  FD  RMAT (13X.2A1.6X.A6.3X.2A10.A5.16X.2A1.6X.A6.3X.2A10 .A5) 

3200  F3RMATI*  ERROR  IN  INPUT  OATA  CARDS  - NO  SYSTEM  NUMBERS  READ  - RUN* 
A ♦ TERMINATED*) 

3300  F)RMAT  (1X.12A10) 

3*00  F3RMATUX) 

EM  0 


non  nnn  o o o n o o o 


*1 


„ .y  . . v THIS  PACK  IS  PEST  QUALITY  rRirrirtur  y 

\r-  *.«r'-  . HS.VM  00l"t  FUWUiHjO)  FOllDC  „ 

SJBKOUTINE  ECHO ( IE  XCH , ENORUN t IRA ) 

THIS  FOJTINE  IS  USED  IN  CONJUNCTION  WITH  SYSTEM  ROUTINE  RECOVR 
TO  R:  INI  TIALIZE  EXECUTION  WITHOUT  RELOADING  IN  THE  EVENT  OF 
EXfiJTIOK  TIHE  FATAL  ERRORS.  EXECUTION  WILL  BE  RESTARTED  ONLY  TO 
PRIN'OU*  COMMON  CAT  A. 


C3NM3N  /INFO/ 

A & RD(8P),  ERR(fe)  , SYS  (40  , ISYS , EON ( 4460 ) , PAGE  ( 25 ,100  > , VAR  ( 10  0 , 21  * 
P EFF.SIOEt 

C CO L ( 1?) , EQN1 ( 960 ) , INPUTS (30) , OPOR ( 120 , 3) , OPR (1 20 , 34) ,PAKC>0,  J)  , 
0 CHAR,II,IFH,1K,IT0,JFH,  J J, JT 0, LINE , LOPR, NCOL ,NEUN, NEQN1 , 

r NIN.NOPR,  NFAR,NVAR,A,  8.C  ,0,  E,F,G,H,I  , J,  , KK  , L , LL  , M , MM,  N , 

F NV ,P,Q,R,S,T,U, V,W,X,V,Z 

3JHHCN  INDEX  M (16  00  II  , EQUATN  (23(i  0)  , OAT  A <4001  , SIZEN  ,S  IZ  M,  SIZEE , 

A SIZER, KEY 
HTEGER 

A CARO, ERR,  SYS,  EON,  PAGE ,VA R,EFF , SIDE , COL, EQN1, OPOR , OPR, PAR, CHA fi, 

B A, P,C,0,E ,F,G,H,P,0,R,S,T,U, V,W,K, Y,Z 
INTEGER  EUUATN.DAl  A, SIZEH, SIZM, SIZEE, SIZER 
31  MEN  SION  IEXCHU7)  ,IRA(1» 

FIGURE  OUT  ERROR  ACORESS  AND  NUMBER 

ER  0RUN*1. 

I>  GO*  *AND (SHIFT  (IEXCH(l)  ,2k)  ,7777773) 

IF  (IAUOR.NE.  0)  GO  TO  10 
I'RROMANDISHIFT  (IRA  ( I ADOR)  ,12)  ,7»> 

IV  00*  =AN0  (SHIFT  (IRA(IAODR) , 30)  ,7777778) 

II IIII-5NMOOE 
G3  TO  20 

1"  irRROMANOCIEXCH(l)  ,77B) 

II  III  I«5HTYPE 

PRINT  COMMON 

20  WR  I Ti  (6,  1000  ) II II II  , I ERROR  , IAOOR, 

A IS  YS ,C  HAR ,11, IFN,IK,ITO, JFN, JJ, JT 0, LINE, LOPR,NCOL, NEQN, 

B NEON1, NIN, NOPR, NP  Afc , N VAR ,A,B,C,D,E,F,G,H,I,J,K,KK,L,LL,N,MM, 

C N,  NN , P, Q, R, S,T,U,®,N,X,V,Z, SIZE H, SIZM, SIZEE, SIZER, 

0 trF,SIDE,KEY,A,B,C,0,E,F ,G,H,I,J,K,KK,L,LL,N,NN,N,NN,P,Q, 

E R,  S,T,U,V,W,X,Y,Z, 

F CA  RO,ERR,  SYS,  EQN,  ( (VAR<I  ,J)  , J«l,2)  ,1*1,100) , C0L,EQN1, 

G INPUTS,  ((OPOR (I  ,J)  ,J»1,3)  ,1*1,12])  , ( (OPR( I , J)  , J*  1,  34)  , 1*1 , 120)  , 
H < ( PAR( I , J ) , J* 1, 3) ,1*1, SO), DATA 

Cl  LL  EXIT 

10 OF  F)RMAT  <*1*//*  JOB  RECOVERED  FROM  ERROR  *A5,02*  AT  AOORESS  *06/ 

A * SIMPLE  INTEGERS*/2(5X»20(I5,1X)/),5X»11(I5,1X)// 

B • SIMPLE  H0LLERITMS*/5X,  3A11/M5X  ,6021/)  5X,S02l// 

C • CARO*/® (SX, 10A11/I // 

0 * ERR* /5X, 6(15, IX)// 

E • SYS* /6 ( JX, 1 OA 11/ )// 

F • EQN*/M»(5X,100A1/),SX,6QA1// 

G * VAR* / 10 0(®X,2All/)// 
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H • C0L*/5X  , 12  ( 15  , IX ) / / 

X * EQN1*/4A(5X»2Q(I5»1X)/)// 

J • INPUTS* /5X  »2Q(I5»1X)/5X,10(I5,IX)// 

K * 0P0R*/1 20(5  X»  3<I5flX)/)// 

L * OPR*/120(5X, 20(15, 1X>/5X,  14(15, 1X>/)// 

N * PAR*/50  (5X, 3(15, IX)/)// 

N * DATA*/40t5X,10All/)  I 

ESO 


WIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
COPY  PUKW1SHH)  TO  HDC 
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SMS  PAW  IS  BEST  QUALHTP1U.CTTCAHLI 

fKVM  CO)  Y FURMSH&D  TO  UDC  __ 


•w  4 

SUBROUTINE  ERROR  (IERKCR)  , RETURNS  ( Aft  AAAA) 

DIAGNOSTIC  PRINT  CUT  ANO  EQUATION  PROCESSING  TERMINATION  SUBROUTINE. 
IER9DR  SIGNALS  THE  TYPE  OF  ERROR  FOR  PRINTOUT  ERROR  MESSAGE 
SELE'TIDN.  THE  NON-STANCARO  ERROR  RETJRN  IS  ALWAYS  USED  TO 
TERMINATE  PROCESSING  OF  THE  OF  FENOIN3  EQUATION. 

SDMNON  /INFO/ 

A CARO  (80  >,  ERR  (6)  , SYS  (4  0)  , ISYS,EQN<*h60>  , PAGE  (25 , 200  ) , VARdOO  ,2)  , 

B tEE,SIOE, 

C CjLI121,EQNI  (960)  , INPUTS  ( 30  ) , OPOR  ( 1 20  , 31  , OPR  < 1 20 , 34)  , PAR  (5  0,3), 
0 CHAR,II,IEH,IK,ITO,jEM, J J, JT0,LINE ,LOPR, NCOL.NEQN, NEQN1 , 

E N1  N»  NOPR, NPAR ,NVAR,A,B,C ,0, E , F , G, H , I , J, K , KK , L, LL , M , MM, N, 

F NN ,P,Q,R, S , T , U, V,N,X,Y,Z 

OWN  INDEX  M (lb  00  14  , EQUATN  (2300)  ,DATA(400)  , SI  ZEM  ,S  12  M,  SI  ZEE  , 

A SIZER,  KEY 
INTEGER 

A Cl  RO,ERR, SYS, EQN, PAGE, VAR, EEE, SIDE, COL, EQNl.OPOK, OPR, PAR, CHAR, 

B A,  B,C«0»E ,F,G,H,P,Q,R,S, T ,U, V , N , X , Y ,Z 
INTEGER  EQUATN, DATA, SIZEM,SIZM,SIZEE» SIZER 
W?  IT:  (6,2000  ) IERROR 

S3  T3  II, 2, 3, 4, 5, 6, 7, 6,9,10,11, 12, 13, 14, 15, 16, 17,  id, 19, 20, 21, 22), 

A IEFROR 

1 M? ITE  (6,1001) 

G)  TO  500 

2 W ? ITE  (6,1002) 

G)  TO  50C 

3 M? ITE (6,1003) 

G)  TO  500 

4 M? ITE  (6,1004) 

G3  T3  50C 

5 W* ITE  (6,1005 ) 

G3  T3  «50P 

6 M? ITE (6,1006) 

S'1  TO  500 

7 M? ITE  (6, 1C07 ) 

3)  TO  500 

8 R7ITE  <6, IPOS) 

G)  TD  50P 

9 W? ITE (6,1009) 

G)  TO  500 

10  MR ITE  (6,1010 ) 

G)  TO  500 

11  WRITE  (6,1011) 

GD  TO  50P 

12  MR ITE (6,1012) 

G3  TO  500 

13  MR ITE (6,1013) 

G3  TO  500 

14  VI?  IT-  (6,1014) 

GD  TO  500 

IK  V? ITE (6,1015) 

GD  TO  500 
16  MR ITl  (6 , 1 C 16  ) 

GD  TO  50C 
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JBC*  Ottf  Y FURNISHED  TO  UDC  


17  M?  ITE  (6,1017  ) OOrY  FURNISHED  TO DDC 

GO  TO  500 

18  MI  IT:  (6,1018) 

53  TO  500 

19  MI  IT*  (6,1019) 

GD  TD  500 

70  MI ITE (6,1070 ) 

G)  TO  500 

71  MUTE  (6,1071) 

GO  TO  50C 

72  MUTE  (6,1027) 

G3  T3  500 

500  MUTE  (6,3000) 

A IS YS, CHAR, II, IFH, I K , ITO, JFM , J J , J \ O , LINE , LOPR.NCOL , NEQN , 

8 NEQNl,NIN,NOPR,NPAR,NVAR,A,B,C,D,E,F,G,H,I, J, K ,KK, L , LL , M ,MM , 

C N, NN,P,Q,R,S,T,U, V,W,X,Y , Z,SIZEM, SI ZM,SI ZEE, SI ZER, 

O EFF,SIOE,KEY,A,B,C,0,E,F,G,H,I,J,  K,KK,L,LL,M,MN,N,NN,P,Q, 

E R,  S, T, U, V , M, X , Y , Z , 

F C& RD,ERR, SYS, EQN, ( ( VAR<I , J) , J=1 , 2 ) ,1*1,100) ,COL,EQNl, 

G INPUTS,  ((OPOR  (I,  J) , J* 1 ,3 ) , 1*1 , 123 ) , ( (OPR ( I , J) , J=l, 34) ,1=1, 120) , 
H <(FAR(I,J),J=1,3I,I  = 1,50),OATA 

RE  TUI  N AAAAAA 

1001  F3 RMA T (*  TOO  MANY  VARIABLES  IN  EQUATION* ) 

1002  FORMAT!*  EQUATION  TOO  BIG  IN  REOUCE3  FORM  FOR  EQN1* ) 

1003  FORMAT  (•  TOO  MANY  PARENTHESIS  PAIRS* ) 

1004  FORMAT (•  MISSING  LEFT  PARENTHESIS*) 

1005  F) RMAT  (*  TOO  MANY  OPERATORS  IN  EQUATION*) 

1006  FORMAT (•  cOULO  NOT  ASSIGN  ALL  OPERATORS  TO  COLUMNS*) 

1007  FO  FMA  T (*  FORCED  ERROR  OUTPUT  - NO  ERROR*) 

1008  F0 RMA T (*  MISSING  RIGHT  PARENTHESIS*) 

1009  F3  RMAT (•  PACK  OR  UNPACK  CALLEO  WITH  NAA  <=0  OR  >*11*) 

1010  F0RMATC*  EQUATION  Oh  FILE  HAS  TOO  MANY  CHARACTERS  FOR  ARRAY  EQN*) 

1011  F) RH1 T (*  EQUATION  DIAGRAM  HAS  TOO  MANY  ROMS  FOR  PAGE  ARRAY*) 

1012  F0RMAT(*  9A0  VALUE  INPUT  TO  SUBROUTINE  VAL*> 

1013  FORMAT (•  INVALID  EQUATION  SYNTAX  - OPERATOR  INPUT  NOT  FOUND*) 

1014  F) RMA  T (*  TOO  MANY  OPERATOR  INPUTS  FOR  ONE  OPERATOR*) 

1015  FORMAT  (•  MISSING  LEFT  SIDE  VARIABLE*) 

1016  F)FMAT(*  MISSING  EQUAL  SIGN*) 

1017  FORMAT!*  RIGHT  SIDE  OF  EQUATION  IS  MISSING* ) 

1018  c0  RMA  T ( • IMPROPER  USE  OF  PRIME  OPERATOR*) 

1019  FORMAT (*  INVALIO  EQUATION  SYNTAX  - OPERATOR  (ANO,  EXCLUSIVE  OR,* 

A ♦ OR,  TIME  DELAY)  NOT  FOUND*) 

1020  FORMAT**  EQUATION  TOO  BIG  FOR  PAGE  ARRAY*) 

1021  FORMAT  (*  BAD  OATA  FOR  UNPACK*) 

1022  FORMAT!*  PACK  OR  UNPACK  CALLED  WITH  38  OR  CC  <=0*) 

2000  FORMAT (*1*//*  ERROR  NUMBER  *13) 

3000  e0  RMA  T (•  TERMINATE  PROCESSING  FOR  T-IIS  EQUATION*/ 

A • SIMPLE  INTEGERS*/2(5X, 20(15, IX)/) ,5X, 11(15, IX)// 

9 * SIMPLE  H0LLERITHS*/5X, 3All/4(5< ,6021/) 5X, 5021// 

C * CARD*/8 <5X, 10A11/I// 

D * ERR* /5X, 6(15, IX)// 

E * SYS*/4(5X,  10A11/)// 

F * EQN*/ 44 (5X, 100A1/) ,5X, 60A1// 

G ♦ VAR*/10  0(5X,2A11/)// 

H * C0L*/5X, 12(15, IX)// 


(AND,  EXCLUSIVE  OR,* 


36 


5JHIS  PAGE  IS  BEST  QUALtTT  FRACTXCA11I 
FROM  COPY  FURNISHED  10  EDO  


EQN1*748  i5X  ,20(15, IX)/)// 

INPUTS*/5X, 20 (I5,1X)/5X,10  (15, IX)// 
GPOR*/120(5X,  3< 15, IX) / ) // 

OPR*  /12  0 (5  X, 20(15,  IX)  / 5X,  14  (15,  IX)  / ) // 
PAR*/50 (5X, 3(15, IX)/)// 
CATA*/*0(5X,10All/)  I 


ooo  non  ooo  ooo  ooo  non  non  ooonooooo 


S J PROUTINE  GETOA T( K£ YFEQ) 
ROUTINE  TO  GET  RANDOM  ACCESS  OATA 
THE  INLY  INPUT  IS  KEYREQ  WHICH  IS 


THIS  PAGE  IS  BEST  QUALITY  PRACTICAAH 
FROM  COPY  FURNISHED  TO  CDC 

FROM  TAPE 7 

THE  KEY  OF  THE  REQUEST EO  RECORO 


THE  ?EC3kO  IS  OUTPUT  IN  CATA(l)  THROUGH  OATA(SIZER).  IF  NO  RECORO 
CAN  3 E FOUNO,  THEN  DATA  IS  SET  TO  QUESTION  HARKS. 


C) MHDN  INOEX  H (16  00  1) , EQU ATN (230  0 ) , DAT  A (400 ) ,SIZEH,SIZH, SI  ZEE , 
A SIZER, KEY 

I H T EGER  EQUATN,OAT  A,S1ZEM,SIZM,SIZE*  , SIZER 

TEST  FOR  REQUEST  TO  OPEN  TAPE7 

Ir (KF  YREQ  »NE  .MHOPEN ) GO  TO  5 

OPEN  T APE 7 ANO  INITIALIZE  BLANK  SOHHON 

SI  7 EH  *16001 
St  ZEE -23CC 

CML  0PENHS(7,IN0EXM,SIZEH,1> 

5> LL  REAOMS(  7 ,EQ  UA  TN  , SI  ZEE  , 9HEQUAT I3NS ) 

5UL  1 NIT  (OAT A,M00, 1H  ) 

SI  ZK  (SIZEH-ll/2 
SI  ZER  *0 
<"  Y*1  H 
Ri  TUR  N 

INITIALIZE 

5 CUL  INIT  (OATA,  400, 10H????7 ??7??> 

FIGURE  OUT  RECORO  SIZE  ANO  SET  SIZER 
S I Z E?  = 0 

0: COOE ( 10 , 10  QQ, K EY REQ I I III  II, KKKKKK , J JJ JJ J 
Ir  ( JjJJJJ.NE.2M  ) GO  TO  10 

EITHER  OISCRETE  OR  SERIAL-OIGITA.  OATA  RECORO  REQUESTED 
SIZER-0 

I7  (II  IXII.EQ.1H-  .OR.  IIIII I.EQ«1H>)  SIZER-16 

Ir ( (IIIIII.GE.1HA  .ANO.  IIIIII.LE.1NZ)  .ANO.  KKKKKK. NE. 1H0) 

A SIZER-16 
GJ  T3  M0 

EITHER  EQUATION  OR  LABEL  OR  EQUATION  SIZE  RECORO  REQUESTEO 
10  If  (JJJJJJ.NE.2HEQ)  GO  TO  30 
EQUATION  REQUESTEO 
D£  COO  E ( 10 , 1100,  KEYREQ)  IIIIII 
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• mf  . l * ^ * i ^ _4  « T v , • 1 vl  * k ' * : ( 

'”0D  2C  JJJJJJ«1,SI2EE 

0- CODEC1C,12«0,EQUATN(JJJJJ  J))  KKKK<K  , LLLLCL 
SI  ZER«ULLLL*0 

20  IF  (IIIIII.EQ. KKK KK K ) GO  TO  40 
SI  ZER*0 
Rr  TWN 

EITHER  LABEL  OR  EQUATION  SIZE  RE30R0  REQUESTED 

30  IF < JJJJJJ.EQ.2HEL)  SIZER«21 
IF  ( JJ  JJJJ.EQ.2HS  ) SIZERxSIZEE 
Ic  (SIZER. EQ.  0)  RETURN 

G:  T OATA 

40  C»LL  KEADMS(7, OATA, SIZER, KEYREQ) 

1000  F3  RMfc  T (4X,A1, 2X,A1,A2) 

1100  FD  RHAT (AB I 
1200  FD  RM4  T ( AS , 12 ) 

?M0 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

SUBROUT  INC  ZNIT  ( X,  N,  AVAL  I WOK  OOP*  fUWUSHU)  TODDC  


C ROUTINE  TO  INITIALIZE  A SELECTED  STORAGE  AREA  (FROM  X<1)  TO  XINM  TO 
C THE  / ALUE  OF  XVAL. 

C 

OT HENSION  X(l) 

INTEGER  Xf XVAL 
03  IP  I»i,N 
10  X<  I ) = X VAL 
TURN 


•'  i'r'yr  Mr.  • ?c  - • "•  ■ 

•..f  * ■ s 's.  t*  m'* 


u«i  ^uAjjin  r 

raojl  COPY  FUivNIijRiuU  TO  DDC 


S J 6R0UTINE  L ABL(  ALPNUM,  VAL, CHAR) , RET URNS (AAA AAA) 

C 

C THIS  SUBROUTINE  MILL  OUTPUT  A TNO  CHARACTER  LABEL.  ALPNUH  SIGNALS 
C NHETiER  AN  ALPHABETIC  (ALPNUH*1)  OR  A NUMERIC  (ALPNUH»2>  IS  OESIREO. 
C VAL  IS  THE  INPUT  VALUE  TO  BE  USEO  TO  GET  THE  LABEL.  CHAR  IS  THE 
C OUTPJT  TNO  CHARACTERS  - LEFT  JUSTIFIED  MITH  BLANK  FILL. 

C 

INTEGER  ALPNUH,  VAL  » RCHAR ,CH AR, R 
I-  (AL  PNUM.EQ.2I  GO  TO  10 

c 

C HAKE  AN  ALPHABETIC  LABEL 

C ( Af  Bf  Cf  0»...|  Z,AA,AB,*..,AZ,BA,8B,...,ZX,ZY,ZZ) 

C I’S  ANO  0*S  ARE  NOT  USEO 

C 

IF  ( VA  L .LT  .1  .OR.  VAL.GT.600)  CALL  ERROR(12)  , RETURNS(20) 

l:har=l*imval-i)/24 

R:HAR«R«VAL-24*(L-1) 

IF (R. GE.9)  RCHAR>RCHAR>1 
Ir  (R.  GE.14)  RCHAR=RCHAR»1 
IF  (L.EQ.l)  LCHAR=55B 

l~  <L.  GE. 2 .ANO.  L.LE.9)  LCH AR=LCHAR- 1 
I-  CL.  GE.15)  LCHARslCHAR^I 

CHAR=  OR (SHIFT (LCHAR,  54) , SHIFT  (RCHAR,  *»8> , 555 5 5 5555 55 55 55 5B) 

RE  TUR  N 
C 

C HAKE  A NUMERIC  LABEL  <01 , 02, 0 3,  ..  . , 10 ,11 , 12,  . . . , 97  ,96,  99) 

C 

10  IF (VAL.LT.l  .OR.  VAL.GT.99)  CALL  ERRORC12) , RETURNS (20) 

l:har*val/io 

R2HAR  = VAL-10*LCHAR 

AR=0R(SHIFT  (LCHAR*33B,54)  ,SHIFT(R;HAR^338,46)  ,55555555555555558) 
RE  TURN 

20  RETURN  AAAAAA 
ENO 
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SJBR3UTINE  P ACK < AA , N A A , BB ,C C) , RETURNS ( AAAAAA) 

C 

C ROUTINE  EITHER  TO  PACK  CHARACTERS  INTO  ARRAY  PACE  (10  CHARACTERS  PER 
C WORO>  OR  TO  UNPACK  CHARACTERS  FROM  AkRIY  PACE. 

C 

C ARRAY  PACE  IS  USEO  TO  STORE  THE  QIAGKAM  PRIOR  TO  PRINTOUT.  PAGE(I,J> 
C HILL  BE  PRINTEO  ON  THE  JTH  LINE  STARTING  ET  COLUMN  10*<X*1)»1. 

C 

C AA  15  THE  CHARACTER <S>  TC  8E  PACKED  INTO  PAGE  OR  UNPACKED  FROM  PAGE  - 
C LEFT  JUSTIFIED  HITH  BLANK  FILL 
C NAA  IS  THE  NUMBER  OF  CHARACTERS  IN  AA 

C BB  IS  THE  PRINTOUT  COLUMN  WHERE  AA  IS  TO  BE  PACKED  INTO  OR  UNPACKED 
C FROM 

C CC  15  TMt  PRINTOUT  LINE  WHERE  AA  IS  TO  8E  PACKEO  INTO  OR  UNPACKED 
C FROM 
C 

C ENTRY  POINT  PACK  PUTS  NAA  CHARACTERS  FRON  AA  INTO  PAGE  AT  RON  CC  AND 
C COLUMNS  BB  THROUGH  BB*NAA-1. 

C 

C ENTRY  POINT  UNPACK  GETS  NAA  CHARACTERS  FOR  AA  FROH  PAGE  AT  ROW  CC  ANO 
C COLUMNS  bd  THROUGH  8B*NAA-1. 

C 

C 

CD  HMD  N /INFO/ 

A CARD  (8  0),  ERR  (6)  , SYS  (40)  , ISYS,  EQN  (4460)  , PAGE  (25  ,200)  , VAR  (100, 2) , 
B EFF,  SIDE, 

C CO  L(  12 ) , EQN1  ( 960 ) , INPUTS  (30)  ,OPOR(  120,3)  , OPR (120 , 34)  ,PAR(50,3> , 
O CHAR,Il,IFM,IK,ITO,JFM, J J, JTO.LINE ,LOPR, NCOL.NEQN, NESNi , 

E NI  N,  NOPR,  NFAR,NVAR,A,B,C,0,E,F,G,H,I,J,K,KK,L,LL,M,HM,N, 

F NN  ,P,Q,R,S,T  ,U,V,W,X,Y,Z 

CD MHON  INDEXN (1600 1) , EQU ATN  (230 0 ) , OAT A (400 ) ,SIZEH,SIZH, SI2EE, 

A SIZER,  KEY 
INTEGER 

A CARO, ERR, SYS , EQN, PAGE ,VA R, EFF, SIDE , COL, EQN1, UPOR, OPR, PAR, CHAR, 

B A,8,C,0,E,F,G,H,P,0,R,S,T,U,V,W,X,Y,Z 
IN  T EGER  EQUA TN, DAT  A, SIZEN,SIZH, SIZES,  SIZER 
INTEGER  AA,BB,CC 

FIGURE  OUT  WHERE  TO  START  IN  PAG-  - AT  PAGE (LL  *KK)  MH  IS  THE 
CHARACTER  POSITION  WITHIN  PAGE(Ll,KK)  TO  START  AT  - HM*1  IS 
THE  LEFTMOST  CHARACTER  ANO  MM«10  IS  THE  RIGHTMOST  CHARACTER 

iriBB.LF.O  .OR.  CC.LC.O)  CALL  ERROR(22),  RETURNS(999) 
L.«(BB-1)/10»1 
MM*MO  C(B9 ,10 ) 

IF  (MM. EQ.  0)  NM«1  Q 

K<«CC 

IF (NAA.LE.O  .OR.  NAA.GE.ll)  CALL  ERROR (9) , R£TURNS(999) 

Ie(LL.GT.25  .OR.  KK.GT.200)  CALL  lRROR(20)  , RETURNS  (999) 

PACK  AA  INTO  PAGE 

THE  PACKING  OF  AA  INTO  PAGE  IS  ACCOMPLISHED  BASED  ON  ONE  OF 
FIVE  POSSIBLE  CONDITIONS  - 
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STATEMENT  CONDITION 
LABEL 

10  H M*  1 AND  MM*NAA<11 

LEFTMOST  NAA  CHARACTERS  OF  AA  GO  INTO  PAGE (LL , KK) 
STARTING  AT  THE  LEFTMOST  POSITION  AND  PAGE (LL*1 , KK) 
IS  N0T  AFFECTEO 
?C  MM=  1 AND  HMfrNAA  = U 

ALL  OF  AA  GOES  INTO  PAGE  ILL , KK)  STARTING  AT  THE 
LEFTMOST  POSITION  AND  PAGE ( LL*1 , KK ) IS  NOT  AFFECTEO 
3*}  MM>  1 AND  MM* NAA <11 

LEFTMOST  NAA  CHARACTERS  OF  AA  GO  INTO  PAGE ( KK, LL ) 
STARTING  AT  THE  MMTH  POSITION  3UT  00  NOT  FILL  OUT 
THE  RIGHT  ENO  OF  PAGE (LL,KK)  AND  PAGE ( LL*1 , KK»  IS 
NOT  AFFECTEO 

40  MM>1  AND  MM*  NAA  = li 

LEFTMOST  NAA  CHARACTERS  OF  AA  GO  INTO  PAGE (LL , KK ) 
STARTING  AT  THE  MMTH  POSITION  AND  00  FILL  OUT 
PAGEILL*  KK)  ANO  PAGE ( LL  *1 ,KK)  IS  NOT  AFFECTEO 
50  MN>1  ANO  MM*NAA>11 

LEFTMOST  NAA  CHARACTERS  GO  INTO  PAGE (L  L » KK) 

STARTING  AT  THE  MMTH  POSITION  ANO  RUN  INTO 
P AGE ( LL* 1 ,KK ) STARTING  AT  THE  LEFTMOST  POSITION 

OE  T ERMINE  PROPER  CONDITION 


1=- 

(MY  .EQ.l 

.ANO. 

MM*NAA.LT.ll) 

GO 

TD 

10 

Ir 

(MM.EQ.l 

.AND. 

MMfrNAA.EQ.  11) 

GO 

TD 

20 

Ic 

(NY  .GT.l 

. ANO. 

MMfrNAA.LT.il) 

GO 

TD 

30 

I- 

(MY.GT.l 

.AND. 

MMfrNAA . EQ. 11) 

GO 

TD 

40 

GD 

TD  50 

10  Y- 1C- NAA 

£ Y COr  E ( 1 0 , 1 1 0 0 , Z ) NAA , Y 

EM  COO  EC  10,  Z,PAGE  (LL,  KK)  ) AA  ,PAGE  ( LL , KK) 

R : T LR  N 

20  °\  GE( LL, KK)  = A A 
R:  TU?  N 

30  Y=MM-1 

NY  = 11  -NAA-MM 

E Y COOE  (IP, 1200, Z)  Y ,NAA  , NN 

EY  root  ( 10,  Z,  PAGE  (LL  , KK) ) PAGE(LL,KK)  , A A,  PAGE  (LL  ,KK) 

RE  TlR  N 

4 0 Y=MM-  1 

£YCODE(10,130C,Z)  Y, NAA 
E Y CODE (10, Z, PAGE  (LL  , KK) ) FA  GE  (LL , KK)  , AA 
R:  TURN 

50  IF (LL*1.GT.25)  CALL  ERROR(20>,  RETURNS(999> 

Y*  MM-  1 

NY*  21 -NAA -MM 

EYCOOE(10,1200,Z>  V , NAA, NN 

EY CODE (20 , Z, PAGE (LL,KK) ) PA  GE (LL , KK) , A A, PAGE (LL  *1 ,KK) 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
Rr  T(R  N row*  cat*  FURNISHO)  TO  LDC  - 

ENTRY  UNPACK 

KILL  START  AT  PAGE(LL,KK>,  MM  IS  THE  CHARACTER  POSITION  WITHIN 
PAGE (LL*KK)  TO  START  AT  - MM*1  15  THE  LEFTMOST  CHARACTER  AND 
MM  = 10  IS  THE  RIGHTMOST  CHARACTER 

E*  TRY  UNPACK 
L.  s (3  B-l)  /10 ♦ 1 
Ml sMOD (B0f 10 ) 

IMHN.EQ.C)  NM*1 0 
K<  »C' 

IF(NAA.LE.O  .OR.  NAA.GE.ll)  CALL  ERR0R(9>,  RETURNS(999) 

Ir  (Ll . GT .25  .OR.  KK.GT.200)  CALL  ERR0R(21) , RETURNS (999) 

Ic  (M,-*NAA.GT.ll  .AND.  LL*1.GT.25)  CALL  ERR0R(21),  RETURNS(999) 

SET  AA  FROM  PAGE 


Y=MM-  1 

EMCOOE(lOilOOOtZ)  Y»NAA 
IE COD  E ( 20  y Zt  PAGE (LL  » KK) ) AA 
R:  TURN 

999  R" TURN  AAAAAA 
1000  FJ  RMA  T ( *(•  II  * X, A#  12  *>♦  ) 

1100  PD  RMA  T ( * ( A*  11  *,  R*  II  *)*  ) 

1200  F) FMAT ( * (A*  II  *,A*  II  *,R*  II  *)*  ) 

1300  FD  RMAT ( * (A*  II  *,A*  11  *)*  ) 

EYO 


ooo  oooo  ooooooooooo  o o o o 
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SUBROUTINE  PARSE,  RETURNS(AAAAAA) 

SUBROUTINE  TO  SEMANTICALLY  ANALYZE  (PARSE)  EQNl(II)  THROUGH  EQN1 ( J J) 


COMMON  /INFO/ 

A CA  RD(0  0),ERR(6)  ,SYS (40) , ISYS,EQN (4460) , PAGE (25 , 200 > , VAR (10 0 , 2) , 
B Es  F, SIDE, 

C CO LC12),E ONI (960) , INPUTS (30) ,OPOR( 120,3) , OPR (120,34) , PAR (50, 3) , 
O CHAR,II,IFM,IK,ITO,JFM,JJ,JTO,LINE,LOPR,NCOL,NEQN,NEQN1, 

E NIN, NOPR, NPAR,N VAR  ,A ,B,C,D,E,F,G,H,I,J,K,KK,L,LL,M,HM,N, 

F NN,P»Q»R,S,T,U, V,W,X,Y,Z 

COMMON  INOEXM(lbOOl)  , EQUATN  (2360)  ,OATA(400)  ,SIZEM,SIZM,SIZEE, 

A SIZER, KEY 
INTEGER 

A CA RO, ERR, SYS, EQN, PAGE, VA R,EFF,SI3E, COL, EQNl,OPOR, OPR, PAR, CHAR, 

B A,B,C,D,E,F,G,H,P,Q,R,S,T,U,V,M,X,Y,Z 
IV  T EGER  EQUATN, DAT  A,SIZ£M,SI7M,SIZEE,SIZER 

THE  BASIC  PROCEDURE  IS  TC  SCAN  EQNl(II)  THROUGH  EQNi ( JJl 
LFFT  TO  RIGHT  - FIRST  FOR  TIME  DELAYS  (TYPE  1 THEN  TYPE  2), 
THEN  FOR  ANO  OPERATORS,  THEN  FOR  EXCLUSIVE  OR 
OPERATORS,  ANO  THEN  FOR  OR  OPERATORS.  AS  THESE  ARE  FOUNO, 

THEY  AND  THEIR  INPUTS  ARE  REMOVE3  FROM  EQNI  ANO  REPLACED  WITH 
1C  CO  ♦NOPR  AND  9999'S.  THE  INFORMATION  FROM  EQNI  IS  PLACEO 
IV  OPR 

INITIALIZE  00  LOOP  FOR  G = 2010,2011,2003,2004,2065 

01  13  C Q=  1, 5 
G=  Q*J  000 

IF  (G.  LE.2002*  G*G»9 
S = II-l 

LDPR=NlN=FxCHAR=0 
C\ LL  INIT( INPUTS, 30, 0) 

LOOK  FOR  OPERATOR  INPUT  (EITHER  VARIABLE  OR  OUTPUT  OF 
ANOTHER,  PREVIOUSLY  PARSED,  OPERATOR) 

20  E*  EU 

CHARxfcQNl  (E) 

Ir  ( CHAR. EQ. 9999)  GO  TO  120 

I; (CHAR.LT.l  .OR.  CHAR. GT. 2 000)  CAL.  ERR0R(13),  RETURNS ( 1000) 
NTN*NIN»1 

IF  (M  N.GT.30)  CALL  EAR0RU4),  RETURNS  ( 1000) 

IF(NIN.EQ.l)  F>E 
INPUT  S( NIN) «CMAR 

LOOK  FOR  OPERATOR  - G IS  CURRENT  OESIREO  OPERATOR 
30  E*EU 

IF  (£.  LE.  J J)  GO  TO  40 
Ir  (LO  FR.EQ.GI  70,130 
40  CH  ARs  EQNI  (E) 

Ir (CHAR.EQ.9999)  GO  TC  30 
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Ir  (CHAR.NE. 20171  GO  TO  50 
C 

C PRIME  FOUND  - RESET  INPUTS(NIN)  »NO  GET  NEXT  CHARACTER 

C 

Ir  (INPUTS  (NIN).EQ.  0)  CALL  ERROR(IS),  RETURNS ( 1060 ) 

INPUTS(NIN)  = -INPUTS(NIN) 

GD  TQ  30 
C 

C TEST  CHAR  FOR  2003,  2004,  2005 

C 

50  Tc  (CHAR.NE. 2003  .ANO.  CHAR. NE. 200*  .AND.  CHAR. NE. 2005  .AND. 

A CHAR. NE. 2010  .AND.  CHAR.  NE. 2011)  CALL  EKR0R(19),  RETURNS  ( 1 00  0) 
C 

C TEST  FOR  DESIRED  OPERATOR  - IF  FOUND  THEN  SET  LOPR  AND  LOOK 

c f:r  hore  operator  inputs 

c 

Ir (CHAR.NE.G ) GO  TO  60 
L ■*  F R=  CHAR 
GT  TO  12C 
C 

C DESIRED  OPERATOR  NCT  FOUND  - IF  PREVIOUS  OPERATOR  (LOPR) 

C WAS  NOT  THE  DESIRED  OPERATOR,  THEN  CLEANOUT  NIN  ANO  INPUTS 

C AnO  START  OVER  LOOKING  FOR  OPERANDS 

C 

60  Ir (LDFR.NE.G)  GO  TO  100 
C 

C LC  PR  WAS  THE  DESIREO  OPERATOR  - CURRENT  OPERATOR  IS  NOT  NOW 

C OESIREO,  THEREFORE  EQN1(F)  THROUGH  EQNi(E-l)  CAN  BE  PUT  INTO 

C O3 R ANO  REMOVED  FRCH  EQN 1 

C 

70  LDPR=t 

ND  PR=  NOPR  ♦ 1 

1-  (NO PR. GT. 120)  CALL  £RROR(5),  RETURNS (1000) 

0*  R (N  OPR,  1)  = G 
0=R  (N  OPR,  2)  = NIN 
DT  8C  P=1 , NI N 
60  03R(NOPR,P*4)=INPUTS(F) 

01  LL  INIT(EQN1(F*1),E“F-1,9999) 

E) Nl( F) =1000  ♦NOP  R 
C 

C CLEANOUT  NIN  ANO  INPUTS 

C 

100  NT  N=F  = 0 

CUL  INIT  (INPUTS  ,30,0) 

C 

C TEST  E TO  SEE  IF  CCNE  WITH  THIS  »ASS  THROUGH  EQN1 

C 

120  IF (E. LT.JJ)  GO  TO  20 
C 

C DONE  WITH  THIS  OPERATOR 

C 

130  CONTINUE 
C 

C EITHER  RETURN  OR  TAKE  CARE  OF  EQUAL  OPERATOR 

C 
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r>  o o o o n 


TAKE  CARE  OF  EQUAL  OPERATOR 


Tc  (E1NK1I.LT. 1 .OR.  EQN1  <1 1 .GT  . 1 j OP  > CALL  £RR0R(15>» 

A K£  TURNS  <10G0) 

Tr (fclNl(2).NE.2006)  CALL  ERROR(16>,  RETURNS (1000 ) 

IF(N(FR.EQ.O  .AND.  INFUTS(l)  .LQ.O)  3AlL  ERR0R(17>,  RETURNS  < 10  00  ) 
Tc  (INPUTS  (1)  .LT.0)  EQN1  (2)  *-2006 
rTW  N 


EFROR  ENCOUNTEREO  - USE  NON-STANOARD  RETURN 


1000  R*  7 URN  AA  AAA  A 

E*D 


; 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
TOO#  QQPY  FURNISHED  TOLDQ  , - 

SJ  BROUTINE  REA04,  RETURNS!* AAAAA, 3B9BBB) 

C 

C ROUTINE  TO  GET  NEXT  EQUATION  FROM  TAPE 4 - LOAOEO  INTO  EQN  - 1 
C CHARACTER  PER  MORO  LEFT  JUSTIFIED  WITH  BLANK  FILL. 

C 

C STANDARD  RETURN  IS  USEO  HHEN  THE  DESIRED  EQUATION  IS  FOUND  HITH  NO 
C PROBLEMS 
C 

C RETURN  AAAAAA  IS  USED  HHEN  THE  EQUATION  IS  TOO  LARGE  FOR  EQN 
C 

C RETURN  BBBBBB  IS  USEO  HHEN  AN  ENO  OF  FILE  IS  ENCOUNTERED 
C 

c 

CDMHON  /INFO/ 

A Cl  FO(80),ERR<6)  , STS  (401,  ISYS, EQN <4460 > , PAGE ( 25 ,200 > , VAR( 10 0 , 2)  , 
9 EFF»SIOEt 

C CD  L(12>,EQN1(  96  0)  , INPUTS  (30  ,OPOR(120,3)  , OPR  (120,34)  ,PAR<50,3), 
0 CHAR,  1 1,  IFF,  IK,  ITO,  JFN,  J J, JTO ,LINE ,LOPR, NCOL ,NEQN,  NEQN1, 

E NI N,NOFR,NPAR,NYAR,A,B,C,0,E,F,G,H,I,J,K,KK,L,LL,M,NN,N, 

F NN,P,Q,R,S,T ,U,V,H,X,Y,Z 

DDMMDN  INOEXM  (1600  1)  , EQUATN  ( 2300)  ,DATA(400)  ,SIZEM  ,S  I2M , SI  ZEE , 

A SIZER,  KEY 
INTEGER 

A CA RO, E RR, S YS, EQN, PAGE, VAR,EFF, SIDE, COL, EQN1,0P0R, OPR, PAR, CHAR, 

B A,  B,C,0,E,F,G,H,P,Q,R,S,T,U,V,H,X,Y,Z 
INTEGER  EQUATN,  DAT  A, SIZEH,SIZM,SIZEE,  SIZER 
C 

C T; ST  CARD ( 1)  TO  SEE  IF  PROGRAM  SHOULD  STOP,  IF  A CARO  SHOULO 

C BE  READ,  OR  IF  EQN  SHOULD  BE  LOADED 

C 

Iff  (CARO(i).EQ.-l)  RETURN  BBBBBB 
10  Ir  (Cl  KO(l).EQ.lHl)  GO  TO  30 
C 

C READ  A CARO 

C 

20  RE  AO<  4,1000)  CARO 

IF (ECF(4) .NE.O.)  RETURN  BBBBBB 
SD  TO  10 
C 

C SEE  IF  A DESIRED  EQUATION  HAS  BEEN  FOUNO 

C 

3"  ENCODE ( IP, 60 0C,H)  CARC ( 2 ) ,C ARO( 3 ) ,Cl RO (4) , CARD (5) , CAR 0(6) , CARD (7) , 
A Cl  RO<8)  ,CAR0(  9) 

ENCODE(10,20C0,U)  H 
ENC 00 E(10, 40 00, V)  H 
ENCO?E(1C,500P,H)  H 
OD  4?  X=1 ,40 

40  Ir  (H. EQ.SYS(X)  .OR.  U.EQ.SVS(X)  .OR.  V.EQ.SYS(X)  .OR.  H.EQ.SYS(X)) 
A GO  TO  60 
C 

C DO  NOT  HANT  THIS  EQUATION  - TRY  NEXT  ONE 

C 

5D  TO  20 
C 

C HANT  THIS  EQUATION  - SET  EFFECTIYITY  (EFF)  AND  EHUX  SECTION 


*8 


i 


THIS  PAGE  IS  BEST  QUALITY  mcriCAlJJl 
FROM  COPY  FURNISHED  TO  DDC  

C (SIDE)  AND  STORE  EQUATION,  HITHOJT  BLANKS,  IN  EQN 

C 

50  E>ICO:E<2,30CO,EFF)  CAPO(74>  ,CAR0(75) 

I7  CC»  KO(80».EQ.1HL>  SIDE*5H  LEFT 
I7  (CARO(80).  EQ.iHR)  SI0E*5HRIGHT 
H-  2 
HZ QN= 0 

60  Ic  (CARO(H>  .NE.iM  ) GO  TO  70 
Ir  < H,  EQ.6  9)  GO  TO  80 
H*  H*i 
G)  TO  60 

70  I7  ( NE  QN . EQ.4  460)  CALL  ERROR  ( 1C ) , RETURNS(999> 

N7QN=NEQN»i 

El  N CM  EQN)  *CARC(H) 

I7  < H.  EQ.69I  GO  TO  80 
H=  HH 
GO  TO  60 
C 

C E>l  O OF  CARO  - READ  NEXT  ONE  TO  SEE  IF  EQUATION  IS  CONTINUED 

C 

80  RE  AO( 4 , 1000)  CARO 

I7  CEDFC4)  .NE.O.)  GO  TO  IOC 
IF  (CA  ROm.EQ.lHl)  GO  TO  90 
H*  11 
G)  TO  60 

90  I7 (NFQN.NE.O)  RETURN 
H -2 

G)  TO  60 

100  I7  (NEON. EQ. 01  RETURN  EBB8BB 
C\RD<1>=-1 
RETURN 

999  RETURN  AAAAAA 
1000  F) RMAT ( 80 Al) 

2000  FORNAT(AA) 

3000  F) RMA T (2Alt 
4000  F) RMAT { A3*0* ) 

5000  FI  RMAT  (A2*00  *) 

6000  FO RHA T (8A1) 

S'!  D 
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